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FRIS PROJECT OVERVIEW 


The Forest Resource Information System Project (FRIS) is a coopera- 
tive effort between the National Aeronautics and Space Administration 
(NASA) and St. Regis Paper Co, (STR), Purdue University's Laboratory 
for Applications of Remote Sensing (LARS), under contract to NASA, will 
supply technical support to the project. 

FRIS is an Application Pilot Test (APT) Project funded by NASA. 

The project is interdisciplinary in nature involving expertise from both 
the public and private sectors. FRIS also represents the first APT to 
involve a large broad base forest industry (STR) in a cooperative with 
the government and the academic communities. 

Purpose 

The goal of FRIS is to demonstrate the feasibility of using computer- 
aided analysis techniques applied of Landsat Multispectral Scanner Data 
to broaden and improve the existing STR forest data base, thereby 
creating the foundation of a dynamic information system. The successful 
demonstration of this technology during the first half of the project 
will lead to the establishment by STR of an independently controlled 
operational forest resource information system in which Landsat data is 
expected to make a significant contribution. FRIS can be viewed by the 
user community as a model of NASA's involvement in practical application 
and effective use of space technology. Additionally, FRIS will serve to 
demonstrate the capability of Landsat MSS data and machine-assisted 
analysis technology to private industry by: 
o Determining economic potentials, 

0 Providing visibility and documentation, and 
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o The ability to provide timely information 
and thus serve management needs. 

The ultimate long terra successful ness of FRIS can be measured through 
future development of remote sensing technology within the forest 
products industry. 

Scope 

FRIS is funded as a modular or Phase project with an anticipated 
duration of three years. The original project concepts were developed 
in 1973, and a formal project plan was submitted to NASA by STR in 
1976. The project officially began in October 1977 after the signing 
of a cooperative agreement between NASA and STR; and after the compel - 
tion of contractual arrangements with Purdue University, 

Organization 

The organization of FRIS is depicted in the chart that follows. 
Since FRIS is a cooperative involving three independent agencies, a 
steering coiranittee consisting of a project manager from each institution 
was formed to provide for overall guidance and coordination. Opera- 
tionally, both STR and LARS have project managers and project staff to 
insure for the timely completion of activities within the project. 

The NASA technical coordinator monitors project activities and provides 
a liaison between the STR and LARS staffs. The solid lines on the 
chart indicate the flow of management responsibility. The dash lines 
reflect the technical and scientific inter-changes between operating 
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1.0 INTRODUCTION 

The material which appears in this report is a reflection of the 
FRIS Project Staff activities for the period 1 January 1980 to 
31 March 1980. This time frame enctmtpasses the fourth quarterly re** 
porting period for Phase III of the Forest Resource Information System 
(PRIS) Applications Pilot Test (APT). Phase III or the System Transfer 
Phase of FRIS is directed at meeting the overall Project goal: 

To document and transfer remote sensing technology developed 
throughout the project that will provide St. Regis with an 
independent operational system^ having Landsat data as a 
significant and viable contributor. 

The most significant ProJect*-wide advance came in late January. At this 
time the St. Regis Corporate management announced that they had accepted 
the FRIS concept, and planned to make available financial resources for 
its support. Corporate acceptance authorized the necessary acquisition 
of equipment, space, and personnel within the Southern Timber lands 
Division to make FRIS a viable entity. Ilils milestone also marked that 
Juncture in the APT where St. Regis personnel assumed the greater bur- 
den of the lead activity and LARS staff assumed a consulting role. 

The unavoidable nine-month delay in the announcement of this 
positive decision may impact the timely completion of the APT. However, 
delay in Project completion is not anticipated to be greater than three 
months, since most System Transfer tasks associated with Phase III were 
undertaken with the premise that Corporate approval of FRIS was forth- 
coming. 

Noteworthy project accomplishments for this last quarterly re- 
porting period Include: 

0 Successful testing of the Landsat 3 reformatting programs at NCC. 
0 Transmission of assembler routines to improve the Landsat 3 
program operating efficiency. 
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0 Compl«tlon of progranaing for the SMOOTH and CHAMCBDBTECTION 
subroutines. 

0 Semi-annual project review at Jacksonville. 

The remaining sections of this report describe these activities in 
more complete detail, /^pendix B contains update timeline charts for 
all tasks. 

2.0 TASK AREA ACTIVITIES 

2.1 Technology Transfer Task 

Technology transfer activities during the past quarter have been 
directed at system installation and operations. The only formal tecli- 
nology transfer activity scheduled during this period was a photo- 
interpretation workshop at Jacksonville. The workshop has been post- 
poned because of time constraints imposed on the St. Regis staff 
associated with the acquisition of FRIS hardware. 

With the transfer of two LARS staff to St. Regis during this 
quarter a smooth transfer of the image processing technology is insured. 
Although this could be a radical departure from traditional technology 
transfer activities, it does underscore the committment of the user to 
remote sensing. Knowledgable individuals within the organization will 
be more successful in disseminating the remote sensing technology than 
would an outsider using more classical approaches. Formal classroom 
sessions, workshops and consulting services are ineffective in answering 
day-to-day operational problems. We are confident that the new FRIS 
employees will be the cornerstone of remote sensing technology within 
St. Regis. 

2.2.1 General System Transfer 

The ultimate operational version of FRIS will be a relatively com- 
plex system composed of three unique subsystems. The three subsystems 
are: 

1) A tabular data base which contains extensive information on 
forest stand conditions. 
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2) A Geographic Inforfflation Systems (CIS) which contains carto- 
graphic* cadastral, and other collateral information* and 

3) An image data base which will be able to both process and dis- 
play information collected by aerospace remote sensors. 

PRIS will be unique in that It will combine digital inforaation from 
all three subsystems in a geographically referenced data base system. 

Tlie system will allow resource managers to interactively monitor and up- 
date conditions of their land base. 

The ability to develop, and make a system like PRIS work is in- 
herent in the level of developiwnt of the three subsystems. The tabular 
data base has been in existence and use for more than ten years. This 
data is resident at the corporate computer facility and accessible via 
remote terminal at Jacksonville. 

The geographic Information system consis'’8 of vendor provided soft- 
ware and hardware. The CIS will be resident on a POP 11/70 minicomputer 
located at the Jacksor.ville PRIS center. In addition to its use for 
creating digitized maps and collateral data, it will be the chief means 
of integrating all forms of digital data. Communications between the 
mini and the mainframe will provide resource managers with the capability 
to access the tabular and image data bases from the mainframe. These 
data will be transferred by land line to the mini where they will be 
combined with graphic map data. Interaction and analysis will occur on 
the mini which will have the capability to produce high quality map out- 
put. 

The core of the image data base will be the LARSYS image processing 
software developed at Purdue. Although this software was developed to 
support research and development programs, it was effectively proven 
during the Demonstration Phase to meet the PRIS objectives. The LARS 
staff objective during this Phase of PRIS is to assist St. Regis staff 
implonent the software at the corporate computer facility. The remain- 
lug subsections of this report addresses the status of this activity. 
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2.2.2 LARSY8 Preprocessing Softwsre Developnent 

LARSYS preprocessing software development Is a task which includes 
three major piecea of software and appropriate documentation. The 
three processors convert digital Landsat data to LARSYS foraat, perform 
systematic geometric corrections of Landsat data and register two images 
of Landsat data. Software documentation will include listing documen* 
tatlon* traditional program abstracts and a user’s guide. 

The first preprocessing system of programs converts digitrl Landsat 
data to LARSYS format. The process used to complete this program was 
to develop; a) functional specltlcation, b) design specification, c) 
implementation plans, and then execution of plans and testing. Functional 
specification refer to program expectations. Design specifications 
refer to program execution. The implementation plan documents the who, 
when and what that relate to tasks required to accomplish the programming 
of the software. Finally, the actual implementation plan that will be 
followed, work to be done, and how the results will be tested. Documen- 
tation in effect takes place throughout this process. The Landsat pro- 
cessor is complete in that it includes documentation, Vrith the exception 
of the users guide. 

The geometric correction is the second major system of pgorams. 

Aside from its multiple corrections of Landsat I and II data, this system 
will be used to rotate the new "P" forr^tted data to true north. Both 
functional and design specification have been completed. Implementation 
will begin during May, 1980. Completion will include most documentation 
and testing and is expected to be wrapped up during July, 1980. 

The last major processor is vhe image registration system. The 
primary purpose of this system is to register two coincident digital 
images as two Landsat digital Inugc data sets. The secondary purpose is 
to provide for the registration <>f any known two dimensional grid to 
another known or defined two dimensional grid. The status of this soft- 
ware is that functional specifications have been completed. Functional 
specifications that define what the image registration system will be 
able to do may be found in Appendix A-1. Tlie overall goal uf this 
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activity is to produce a maintainable, efficl«>ntt system which is 
modulised, well documented, and easy to use. Both cubic and nearest 
neighbor interpolation will be available* Locations may be approximated 
by up to a third order bicubic polynomial. 

Design specifications and implctocntation plans are to be completed 
by mid-June, l'*80. Some implement .it ion has already begun. Cubic inter- 
polation will follow the algorithm described in Appendix A-2. Control 
cards have been carefully cc.tstructed to cover all Tunctional require- 
ments as well as simplicity of use. Other main image registration 
section elements and design specifications are still being finalized. 

The multifit lesst squares analysis section has been initiated. The 
cross-correlation section will be started in June, 1980. 

Documentation is the last major effort of the LARSYS preprocessing 
software implementation. The three types of documentation will be pro- 
gram listing documentation, program module abstracts, and user documen- 
tation. The first type, listing documentation, is most important to 
system analyst personnel who will maintain these programs. Because of th 
this key concern for maintainability, a draft of a standard for listing 
documentation has been generated for this project. Features are that the 
leading comments in a listing will contain such key information what the 
program does, what the inputs and outputs are, and all global and local 
variable descriptions. Legibility of code is specifically emphasized 
(refer to Appendix A-3). 

The total 1.ARSYS preprocessing effort is progressing at a steady 
pace. The Landsat reformatting is virtually complete while geometric 
correction is next closest to completion. Image registration will have 
the most effort applied to it during the last two quarters of the project. 

Another important preprocessing transfer activity involves the 
future potential use of fully corrected, P-format, data a*'ailable from 
EDC. The availability of P-type data to FRIS will eliminate mud. of the 
front-end preprocessing currently required prior to image classification. 

The discussion that follows gives preliminary results on the use of fully 
corrected Landsat 3 data from the Picayune test sixe in southern Mississippi. 
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The fully geonetrical ly corrected Landsat MSS frames acquired for 
the forest resource data base are placed In a specific projection and 
orientation. This makes possible a one-to-one correspondence between 
earth coordinates and row column pixel locations in the data. Having 
such a relationship for each frame will enable resource polygons on 
maps to be automatically related to row column locations in the data. 
Visual searching in the imagery would then be unnecessary once comer 
latitude, longitude, or UTH coordinates were known. A program is being 
developed to enable user conversion of coordinates and some of the de- 
tails are included here. 

The fully corrected MSS data are placed in a Hotlne Oblique Mercator 
(HOM) projection and in the future they will be placed in the Space 
Oblique Mercator (SOM) projection. These projections are discussed in 
Appendix D of the new Lcndsat User's Handbook. The scale distortions of 
these projections is very small (1:10,000); thus a linear transformation 
can accurately be used to relate points In the frame. The earth is 
divided into zones of latitude and within each zone the corrected frames 
have a constant azimuth. The zone covering the areas of Interest here 
is zone 2 with latitude range 23*^ N to 48** N and the zone azimuth is 
14.3394993®. The pixel scale of the fully corrected data is 57 meters in 
both directions. 

The software in its present form utilizes a latitude-longitude to 
Universal Transverse Mercator conversion program to transform user input 
latitude- longitude coordinates first to UTM. Then a linear conversion 
is made to line column using the expressions: 

LINE = CLINE + DLINE 

COL « CCOL + DCOL 

DLINE = (-DELEAS . SIN (ALPHA) - DELNOR . COS (ALPHA) )/ 57. 

DCOL => (DELEAS . COS (ALPHA) - DELNOR . SIN(ALPHA)//57. 

DELEAS * EAST - CEAST 

DELNOR * NOR - CNOR 
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where: CLINE» CCOL are the center line and column of the frame. 

CEAST, CNOR is the UTM easting and northing of the center 
point. 

EAST, NOR are the UTM easting and northing of the point to 
be transformed to LINE* COLUMN. 

The conversion program (L.OCPNT) is being developed for interactive 
terminal use and will require typing In the frame center latitude and 
longitude; then the user enters any number of latitude-longitude points 
In the frame he wants to convert. Problems are currently being encountered 
In testing the program on the Picayune frame with Inaccurate results. 

Four test points were taken from the Nicholson and Dead Tiger Creek 
quadrangles in the Picayune frame and the latitude-longitude coordinates 
were input and the output line and column were observed. The input 
parameters are a part of the problem. A latitude and longitude are 
given as the frame format center; however, it was uncertain what exact 
line-and-column number corresponded to this. The bias observed at one 
of the test points was removed and the resulting center line column was 
taken as the format center. Thus, there is no error at this point. At 
the other three points, errors were observed. The results arc listed 
in Table 2.2.2.1 . 

TABLE 2. 2. 2.1 Coordinate Conversation Tests for Picayune Frame. 

Format Center: 30.18° N. , 89.52° W. Center Line, 

Col: 1518,1796. 


Lat. 

Test 

Long. 

Point 

Line 

Col. 

Estimated 

Line 

Point 

Col. 

Error 

Line 

Col. 

30.375 

89.625 

1189 

1518 

1189 

1518 

0 

0 

30.5 

89.75 

987 

1265 

999 

1285 

12 

-7 

30.375 

89.5 

1150 

1725 

1141 

1724 

-9 

-1 

30.5 

89.625 

948 

1473 

952 

1463 

4 

-10 


Causes for these errors are being Investigated. 
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2.2.3 LARSYS Transfer 

The LARSYS transfer task Involves only those processors associated 
with image classification. The bulk of this software was transferred 
during the previous two quarters. The software transferred contains 
elements of LARSYS ver. 3.1 and LARSYSDV. LARSYSDV includes new pro- 
grams which are under development and not part of the ver. 3.1 documen- 
tation. The programs were transferred in card image format on 9-track 
computer compatible tapes. Copies of tape listings and user documen- 
tation were also provided. 

The tasks facing St. Regis personnel are to convert the programs 
which currently run on an IBM 3031 operating under VM to an IBM 3033 
operating un.^er MVS. That is the LARS computer operates as a virtual 
machine while the St. Regis computer operates as a batch machine. This 
means that the LARSYS programs are not directly compatible between the 
two IBM machines. 

St. Regis staff are required to make certain changes to the LARSYS 
software. A summary of the necessary changes to the software appears 
below: 

A. Add the function COPYTAP, this allows the data to be read from 
tape to disk and stored on disk. St. Regis has a disk based 
system while LARS is tape based. 

B. Replace command language with ROSCOE. Due to the operating 
system differences between the two machines, the command 
language has to be modified. ROSCOE Is a software package 
that permits St. Regis users to Inltiace jobs from remote 
sites. This will replace CMS currently used in LARSYS to 
perform similar functions. 

C. LARSYS contains some bookkeeping routines that will be deleted 
because these functions are already handled by St.rRegis. 

D. All non-standard file handling routines in LARSYS will be re- 
placed to meet St. Regis computer software conventions. 

E. All tape handling routines will be modified to deal with disks. 


F« Machine dependent assembly language routines will be eliminated 
where feasible. 

Implementation tasks facing St. Regis staff Include: 

1. Program compilation from tape. 

2. Creation of disk files. 

3. Modification of software for compatibility to St. Regis 
machine* Including elimination of bookkeeping, assembler 
routines and modification of tape callable routines. 

4. Creation of ROSCOE modules, 

5. Development of links to CIS, and 

6. Develop St. Regis/LARS user documentation. 

LARS staff are available for consultation and debugging as needed. 
Our experience during this past quarter was that very little assistance 
was requested by St. Regis personnel. Implementation of these software 
are progressing with very few problems. This is most likely due to, 

a) the level of documentation provided with the LARSYS software, and 

b) the knowledge of the staff Involved with the implementation. 

2.2.4 Programming Additions 

The LARSYS software packages were originally designed to process 
digital multlspectral scanner data in a research environment. Periodic- 
ally, modifications and embellishments have been added to LARSYS support 
packages to improve interaction with the human component of the analysis 
activity. Since FRIS is a user oriented, operational system there were 
certain additions required to improve user efficiency. There have been 
a number of additions to the LARSYS software since the midpoint of 
Phase II. The two newest additions reported this quarter are significant 
because they directly affect FRIS requirements. The two new program 
additions are SMOOTHRESULTS and CHANGEDETECTION. 

SMOOTHRESULTS is a post classification processor designed to emulate 
the human action of creating a mapping cell. Mapping cells are the 
basic component of timber type or operating area maps. The theory 
behind the mapping cell is simply that areas less than a minimum size. 
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say five acres» are Ignored for map drawing purposes and included as 
part of a larger population. Therefore, a two or three acre inclusion in 
a type would be ignored when the map is created. 

The human quickly handled these small inclusions when making a 
type map. A Landsat classification, however, will display most inclu- 
sions that fall within the scanner resolution. These will result in a 
salt and pepper effect on classification output. A situation that may 
accurately portray the cover composition uut which is often not appealing 
to land managers who are used to working with "clean" (no salt and 
pepper) maps. 

SMOOTHRESULTS allows the analyst to define a mapping unit and produce 
a classification results map which does not exhibit a salt and pepper 
pattern. The processor scans a LARSYS Classification Results File and 
replaces groups of classified points (cells) with the dominant class from 
that group. The analyst has the option to specify the size of the cell 
(CELLSIZE card), class numbers which are to be replaced (PRIORITY card) 
and weighting factors for each class (WEIGHTS card). The output from 
this function is to tape or disk in LARSYS Classification Results File 
format. Figure 2. 2. 4.1 is an example of a classification result which 
shows output both before and after use of the SMOOTHRESULTS processor. 

SMOOTHRESULTS was debugged and tested during this quarter. An 
additional option which allows the analyst to define new classes which 
are mixtures of old classes was developed and Is being tested. The 
control card reference file and program abstracts are included in 
Appendix A-4. 

The other processor that was upgraded for addition to the FRIS 
package of software which is being transferred to St. Regis is 
CHANGEDETECTION. This is another post classification processor, and is 
designed to make comparisons between classification results. This 
processor is intended to be used to compare two anniversary Landsat 
classifications which have similar class structures. The resulting pro- 
duct of this comparison is a LARSYS results tape containing "change" 
classes. 
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Change classes are designated by the analyst and are in the form 
where: 

Pine (time 1) goes to Non-Pine (time 2), and 

Non-stocked (time 1) goes to Stocked (time 2). 

Optimally, Landsat data is of an anniversary nature, that is the data of 
collection for both dates is nearly coincident but chronologically a 
year or more apart in time. Present requirements of the CHANGEDETECTION 
program are that the Landsat scenes be precision registered. Independent 
classifications are generated for time 1 and time 2. The analyst is 
careful to Insure that class structure, that is the various spectral 
groups that comprise the information classes is similar. Once the classi- 
fications have been generated, CHANGEDETECTION is ran, and an output 
similar to figure 2.2.4. 2 is produced. Tabular information which 
Indicates the amount of change in acres percent of area by class can 
also be produced. 

During this past quarter, program abstracts for CHANGEDETECTION 
were completed, and drafts of user documentation were prepared. Copies 
of this material appear in Appendix A-5. 

2.3 Classification Evaltiation Review 

Traditionally, performance of multlspectral scanner data classifi- 
cations have been assessed by the evaluation of test fields. Test fields 
vary in size from single pixels to multiple pixel blocks which are 
located randomly or systematically throughout the classification. The 
number of test fields can be statistically determined so that results 
can be stated for a given confidence range. A priori Information is 
used to help the analyst define the number of sample test fields needed. 

Ideally, test fields should be homogenous, that is they should 
represent a 'pure* cover type. Recent trends in evaluating Landsat 
classifications have been toward selecting single pixels or blocks no 
larger than four pixels. This procedure has apparently worked well when 
agricultural or general land use classifications are being evaluated. In 
fact, one may go so far as to assume that the quality of the ground 



CLASSES 


Any class chanfje 
to nonstocked 


Nonstocked 
Same both dates 


Low densUy pine 
Same bcth dates 


Pine and Slash/Cypress 
Same both dates 


*1ixed Pine/hardwood 
Same both dates 


Nonstocked to Any vegetation 


Figure 2. 2. 4. 2 Output of CHANCEDETECTION program; classes shovm are 
change classes. 
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reference data for these types of classifications has a direct bearing 
on the success of this approach. 

In dealing with wildland type classifications these traditional 
evaluation approaches are not as effective. Minimally, these procedures 
are more difficult to implement and consume more time because Identifying 
a 'pure' cover type pixel becomes more difficult. Even If a sufficient 
number of 'pure' pixels are Identified performance results can be mis- 
leading, unless the analyst can carefully relate the pixel location to 
concurrent aerial photography. Naturally, the analyst must be able to 
Interpret the photography to determine what test pixels are valid. This 
Is a time consuming process and often impractical because concurrent 
photography Is not always available. 

The most readily available form of ground reference information used 
to evaluate natural resource classification performance Is a cover type 
map. Such maps are developed as aids to natural resource managers and 
therefore tend to represent a cover type as a management entity rather 
than as a spectrally homogenous area. This is an important distinction 
because a type map may contain a spectrally heterogenous mix of pixels 
that from a managers viewpoint form a homogenous mapping unit. If the 
analyst assumes that the resource managers definition of homogenous is 
similar to the evaluation criteria of purity the performance results 
will be poor, when in fact performance may be good. It is for precisely 
these reasons that the traditional evaluation methods are not adequate 
for determining classification performance for wildland classification. 

The classification evaluation procedures used for FRIS were of Che 
traditional variety. Based on this experience we recommend implementa- 
tion of procedures Chat operate with the benefit of an Interactive 
Graphics System (IGS). The following chart (Figure 2.3.1) illustrates 
the general procedure used to evaluate a classification and how this 
procedure differs when an Interactive graphics system is available. 

Previous Methodologies 

Various strategies have been employed to sample St. Regis Ownerships 
for classification accuracy evaluations. During the FRIS demonstration a 
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formula was developed which determined the sample size needed for a 
specified confidence range around the classification accuracies to be 
evaluated. Then a systematic or random sampling scheme (with test fields 
of one pixel) was applied to three test sites. The systematic samplli^ 
method was the preferred method for the analyst when applied carefully 
and with full knowledge of its cyclical nature. The systematic sampling 
method was precise and had less human error involved than the random 
sampling approach. 

Various test field sizes have been employed to evaluate St. Regis 
classifications. Initially » test fields of four pixels each in a 
systematic array were used for classification evaluation. Heterogeneous 
test fields were dropped as well as test fields which fell on irr>scon- 
cllable map boundaries. Ultimately only A2.7Z of the original fields 
were utilized to assess classification performance. Future evaluations 
utilized a single rather than a multiple pixel test field, thus increasing 
the number of samples and decreasing Che man-time that was Involved with 
the multiple pixel systematic sample. 

Previously, as illustrated in Figure 2.3.1, without IGS, the 
selection of test fields for evaluation presented many problems, that is: 

1) What confidence was required on the accuracies produced? 

2) How would the test fields be selected (randomly, systematically, 
stratified random samples, etc.)? 

3) How large would each test field be (one pixel or a block of 
pixels)? 

Since the IGS mechanises this step in the procedure, thus speeding up 
the process dramatically, the whole area can now be used for evaluation. 
Hence any decisions involving confidence limits, types of samples, or 
size of test fields are eliminated. 

Another problem does remain, however. The problem involves handling 
border pixels in classification evaluations. 

The rest of the steps Involved in evaluating classifications can 
all be mechanised within the Interactive graphics system uaployed at FRIS. 
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Evaluation 

Procedure without I6S 


Evaluation 
Procedure with I6S 


#conf1dence% 
' required * 


determine 

number of 
test fields 


\ 


> 

select test 
fields 

a) randomly j 

b) systematically 


all fields 
canTe TficTuJed 
automatically 


♦ PRINTRESULTS 
map test fields 


not necessary 


OVERLAY 
test fields 


can be automatically 
done 


RECORD ^ 

Landsat and 

Inventory categories { 

■ automatically 

’ recorded and tallied 

for classification 
j accuracy 

Classification Accuracy / 

and confusion table 


Figure 2.3.1 Comparison of Classification Evaluation Procedures with 
and without an Inter‘<^:tlve Graphics System (ICS). 
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Thus* many man hours and the resulting human errors can be eliminated* 

As stated earlier and Illustrated In Figure 2.3*1 the IGS improves 
classification evaluations In two key areas: 

0 The whole area or a very large sample of the area can be 
utilised for evaluation. 

0 The system can be virtually* completely* atuomatic* eliminating 
many man-hours of work* 

Uith these Improvements come new and In some respects different problems* 

First of all, sir re polygons* not single pixels or blocks of 
pixels* will be used for evaluation* criteria must be developed for 
labeling an heterogeneous polygon. The proposed criteria for FRIS Is: 

0 75Z or more pixels pine -> "pine** polygon 

0 75X or more pixels hardwood » "hardvrood" polygon 

0 otherwise the polygon Is labeled as "other” or "nonstocked" 

These criteria may change as the spectral characteristics of these cover 
types are studied further and/or more classes or subtypes become 
separable* 

Another problem involves the polygons' border pixels and how they 
are to be treated with respect to the classification evaluation. In 
ocher words* If a pixel lies on the boundary between two or more polygons 
(Figure 2*3.2) of what polygon should it be considered a member, or 
sliould It be considered a member of any polygon? There are three 
possible solutions for handling this problem: 

1) Include all border pixels in whichever polygon the IGS overlay 
has assigned them* 

2) Exclude all border pixels from the analyses. 

3) Escclude only some border pixels from the evaluation. 

The first solution* including all border pixels* would give a com- 
plete* unbiased evaluation of the area. In this case, however* any small 
error In registration could substantially effect the evaluation, causing 
border pixels to be assigned to the incorrect polygon. Also, this method 


T 


r 

I 


Ll 

LI 


- interior pixel 

- border pixel 


Figure 2.3.2 Overlayed area illuatratlng border pixela and interior 
plxela. 
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atsuMS tnat apactrally a pixel should fall in the claee which covers 
the aajority of the pixel. In fact» however » this is not necessarily the 
case. This nethod would be extrenely feat and require very little 
programing and hence could serve as a first cut at evaluating a ciasal~ 
ficatlon. 

The second solution* excluding all border pixels, vrnuld give an 
Incomplete, (not all pixels Included), biased evaluation of the area. 

If, however, the assumption can be made that the border pixels of any 
polygon have tlM same distribution of classes as the in^ <hI pixels, 
then this solution becomes \mbiased. Also, this method . ^:cet the 
asK>«mt of reglstrati<’M. and mapping error measured in the evalvacion. 

Some programming, supplementing the IGS will be necessary i-j - 'p lament 
this solution. 

The third solution is essentially the same as the second solution 
except that those pixels which are on a burder between two spectrally 
undifferentlable polygon types would be included in the evaluation. That 
is. If a pixel lies on the boundary between two pine plantations which 
differ by a spectrally unmeasurable criterion (i.e., 2 years age 
difference), that pi cel would be included in the evaluation. This 
solution has the same criticisms as the second solution (excluding all 
the border pixels) except that more of the area is being evaluated thus 
giving a somewhat better total evaluation. With this solution, howev«^r, 
the assumption that the remaining border pixels of any polygon reflect 
the distribution of those included, may not be as reasonable. 

Another problem involved with using an interactive graphic system 
for classification evaluation concerns the size of the polygons evaluated. 
If a small number of pixels (e.g., 1, 2 or 3) are used to assign a poly- 
gon’s class, the assigned class could easily be in error. Thus a 
criterion, such as 10 internal pixels, should be developed, indicating 
the smallest polygon to be evaluated. 

The interactive graphics system which FRIS will employ for classlfi- 
caticvi evaluation will facilitate and improve this procedure greatly. 
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The new procedure, once It is implemented, will become virtually auto- 
matic. Many previous sampling consldersatlons will become unnecessary 
due to the graphics and overlay capabilities of the new system. 

2.4 MANAGEMENT 

The FRIS management activity oversees day-to-day project operations 
and is responsible for all technical and fiscal project reports. Status 
of all major System Transfer Tasks are shown in Exhibits 1 through 4 of 
Appendix B. A brief discussion of these status charts follows. 

Technology Transfer Task 

All computer-aided analysis workshops and short courses have been 
successfully completed. Activity A3, the photo-interpretation short 
course has been extended twice and in all probability will be eliminated 
as a Phase III activity. Time constraints on St. Regis staff involved 
with the Installation of FRIS systems prohibit selection of a time during 
the remainder of Phase III when this course can be given. The possibility 
of scheduling this activity during Phase IV will be Investigated. 

Consultation (activity B) associated with LARS software implementa- 
tion will be continued through the end of Phase III. Onforseen delays 
associated with St. Regis' acquisition and installation of vendor hard- 
ware and software at Jacksonville are the reasons for this shift. 

User documentation, especially activity C2, is approximately 66% 
complete. Part of this activity will carry over to Phase IV. 

Remote Terminal activities will cease on 30 June as projected. 

LARSYS Transfer Task 

The first two activities of this Task, Planning and Transfer are 
complete. Activity C, Consultation and Debugging is extended to 30 June. 
At this time all LARSYS processors should be operational at NCC. 

St. Regis staff will be primarily responsible for debugging activities 
with LARS staff assisting as required. 

r 


I 
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Documentation (activity D) will be extended into Phase IV. 
Supplen^ntal timeline charts to Exhibit 2 in /^pendix B which will 
Identify specific documentation tasks are included. 

Activity E (Test and Evaluation) is anticipated to occur during the 
last quarter of Phase III. St. Regis staff will be primarily responsible 
for this activity. 

Preprocessing Transfer Task 

The Planning and most of the Program Refinement activities (B1 and 
B2) are complete. Activity B3, Image Registration, will carry over to 
Phase IV, as will that portion of activity C, Program Transfer, and 
activity D, Consultation and Debugging. The image registration software 
will be transferred early during Phase IV. 

Documentation, activity E, is well advanced for the reformatting 
and geometric correction software. The image registration documentation 
is being prepared in parallel with the software programming. Because 
of this, documentation will be completed during Phase IV. 

Test and Evaluation, activity F, is dependent on St. Regis staff 
and their use of the systems. The activity will not be entirely com- 
pleted until transfer and implementation of the image registration soft- 
ware. 

Support, activities G, will be completed by the end of Phase III, 
with the possible exception of G3. Reformatting Operations Procedures 
are closely associated with communications capabilities between the 
mainframe and the mini. LARS staff will assist St. Regis personnel in a 
consulting role as needed to support this activity. 

Management Task 

Management activities are up to date as indicated by the status 
chart. Noteworthy among these activities is a new task being developed 
under Information Dissemination. Specifically, we are developing a FRIS 
color brochure for wide dissemination by the three agencies. Content 
and format of this product have been discussed and a cost proposal to 
cover this work is being prepared. 
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impend ix A-1 

Image Registration Functional Specifications 


An image registration capability has been determined to be a 
necessary part of the FRIS III image preprocessing software. Image 
registration is general enough to mean grid to grid transformation. 
Thus, while the system is designed to register two coincident Landsat 
scenes, registration to alternate grid systems may be accomplished 
with this software as well. Functional specifications will be as 
follows: 

I. Purpose 

A. Primary: Registration of two coincident digital images as 

two Landsat digital image data sets. 

B. Secondary: Provide for the registration of any known two- 

dimensional grid to another known or defined two-dimensional 
grid. 

II. Input images are assumed to be in LARSYS format. 

III. Checkpoint Acquisition 

A. Manual checkpoint acquisition is possible. 

B. Cross-correlation of two coincident digital images may be 
accomplished by implementation of a numerical integration 
image correlator. 

C. Control may be by set line and column intervals. 

D. Alternate control will be from a set of inputted control 
correlation point locations where a cross correlation is 
desired, i.e., arbitrary point by point correlation. 

IV. Registration transformation 

A. Coefficient determination will be calculated for affine, 
biquad, and bicubic transformation. 

B. Transformations through bicubic will be implemented for 
the registration transformation. 

C. Block registration technique will be utilized. 

1. Optimum rectangular block size will be determined for 
biquadratic and bicubic registrations. 

D. Interiors of all blocks will be registered with an affine 
or linear transformation. 
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Appendix A-2 

Cubic Interpolation Used in the Image Registration System 


The algorithm used in the current image registration system for cubic 
interpolation of data values is based on a thrid order Lagrange inter- 
polation. ^nie general Langrangian Interpolating polynomial for three 
dimensions is: 
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The image registration system uses the above equations with m ^ 3, n > 3. 
Therefore, we need m-*-l>4 different X^ values and 0+1=4 different Y^ 
values. The X^'s and used are 0,1, 2, 3 and 0,1, 2, 3. Then the general 

eguaticm reduces to: 

P33(X,Y) = LQ(X)LQ(Y)f(0,0)+Lj^(X)LQ(Y)f(l,0) + 

L2(X)LQ(Y)f(2,0)+L3(X)LQ(Y)f(3,0) + 

Lq (X) L 3 (Y) f (0 , 1) +L3 (X) L3 (Y) f ( 1, 1) + 

L2 (X) L^ (Y) f (2 , 1) +L3 (X) L^ (Y) f ( 3 , 1) + 

Lq (X) (Y) f (0, 2) +L3 (X) L2 (Y) f (1, 2) + 

L2(X)L2 (Y)i( 2,2)+L3(X)L2 (Y)f (3,2) + 

Lq (X) Lj (Y) f (0, 3) +L3 (X) L3 (Y) f (1, 3) + 

L2<X)L3(Y)f(2,3)+L3(X)L3{Y)f (3,3) 
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wheret 


IjW - 

(X-1) (X-2) (X-3) 

X - 6X +11X -6 

(0-1) (0-2) (0-3) 

-6 

L, (X) = • 

(X-0) (X-2) (X-3) 

X^ - 5X^ + 6X 

(1-0) (1-2) (1-3) 

2 


L (X) = (X-0) (X-1) (X-3) X^ - 4X^ + 3X 

(2-0) (2-1) (2-3) " — 

r IV) - (X-0) (X-1) (X-2) _ X^ - 3X^ + 2X 

3' ' - (3-0) (3-1) (3-2) " 6 


2 Uid L^(Y)'s have the same equations with Y substituted for X 
and f(X,Y) is the dsta value associated with pixel (X,Y). 


lb save computation time, the L; ' s agre calcualted according to the above 

equati(»is for specific points in the (X,Y) grid. These points were chosen 
at quarter pixel intervals as shown in figure 1. The calucalted Lj|^(X)<s aire; 



L,(X) 

Lj^(X) 

1.2 1« 

(X) 

X » 1.00 

o 

• 

o 

1.0 

0.0 

0.0 

X « 1.25 

-0.0546875 

0.8203125 

0.2734375 

-0.0390625 

X = 1.50 

-0.0625 

0.5625 

0.5625 

-0.0625 

X » 1.75 

-0.0390625 

0.2734375 

0.8203125 

-0,0546875 

X * 2.00 

0.0 

o 

t 

o 

1.0 

0.0 


Ihe same table applies for Y=1.00, 1,25,1.50, 1.75, 2.00. 

In the image registration process, an input point A (see Figure 1) is 
approximated to its nearest quarter pixel, To calculate the data value 
associated with A, the Lagrange polynomial coefficients for that quarter 
pixel location are used in the P^^(X,Y) equation. To further save on 
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computation« the products L. (X)L. (Y) for all combinations of the quarter 
pixel locations ((1.0, 1.0), (1.25,1.0), (1.50,1.0), (1.75,1.0), (2. 0,1.0), 
(1.0,1.25), (1.25,1.25), etc.) have been stored in a table. Ihen when 
P 2 j(X,y) is calculated, a table loo)c«q> locates the appropriate L^(X)L^(X) *s. 

When this algorithm was inplemented for cubic interpolaticxi of data values, 
it was determined that the error introduced by this method of using 
discrete intervals versus continuous intervals was negligible. It was 
negligible because the intervals involved were quarter pixels and the 
final data values %irere integer values between 0 and 255. 


References: 

"Hultiteiqporal Image Registrations of Multispectral LANDSAT Data of 
Finney and Ellis Co.'s, Kansas by Nearest-Neighbor and Third Order 
Lagrangian Intexpolatics) Methods." Prepared by Charles R. Smith, LARS, 
September 20, 1976. 

Source listing of OVERLA subroutine used in current Image Registration 
System. 


Figure 1» 4x4 Data matrix surrounding point 

to be interpolated (point A). Example: 
Since point A is nearest grid coordi** 
nates (l.S, 1.75), the Lagrange coeffi- 
cients for this x and y are taken from 
the table and used in the interpolating 
polynomial. 
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Appendix A-3 

Reformatting Documentation Standards 


Ft'eface 


Ihis guide si^plements the lARSYS Standards Report Section III 
Programming Standards. Programmers writing software for the Reformatting 
group should read the LARSYS report as well as this guide; wherever this 
guide conflicts with the LARSYS report, this guide should be followed. 
Programmers should take particular note of the paragraphs in the LARSYS 
Standards Report Section 111 on Assembler and EXEC organizations and 
coinnents, and on programming techniques. 

The main enphasis of the guide is on the uocumentation of program 
source code. Program logic must flow downward, and comments must 
reflect that flow. Within the source code, all global and local variables 
must be identified in variable description lists. The source code also 
must contain a general description of the algorithm used and input/output 
requirements. Specific coding and commenting practices are recommended 
for inproving the legibility of source code. 

This guide contains the following information. 

I. Documentation Outside of Source Code Listings 

II. Documentation Within Source Code Listings 

A. Overall System Standards 

B. Layout of Individual Routines 

C. Comments Within the Body of Routines 

Appendix A Exanple Control Card Description 

Appendix B Exanple LARS Program Abstract 

i^pendix C Example Software System 


j^pendix D Example Block Data 


DOCUMENTATION OUTSIDE OP SOURCE CODE LISTING 


Any program with a control card reader must have a separate 
description of its control cards. The description must include 
all keywords and all parameters with an indication of which keywords 
and parameters are required and which are optional. All default 
values must be indicated. It is also useful to include one or two 
sanple control card decks. For an example of a control card 
description, see Appendix A. 

B. Any program designed for use by non-reformatting staff should 
have a user's guide. This guide should include several example 
user sessions. 

C. Any program using routines that depend on non-trivial algorithms, 
calculations, or data structures must have an ad>stract. The abstract 
may be for an entire system or for specific subroutines. The 
abstract must describe the edgorithms, calculations, and/or data 
structures in sufficient detail for a person unfamiliar with the 
source code to understand the inplementation. For a major program, 
it may be appropriate to have two levels of documentation 
abstracts. One abstract would be directed at the interested user, 
emd the other at the programmer responsible for progr6un maintenance. 

^ For an exanple of a program abstract, see Appendix B. 

i 

t 
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XI. DOCUMENTATION WITHIN THE SOURCE CODE LISTINGS 
A. Overall System Standards 

1. Bach Routine must flow logically downward. See J^pendix C 
for exaiqples of routines that flow logically downward. 

2. The names of all routines for a specific software system 
must have the same three*letter prefix. The last three 
letters should be unique for each routine and represent the 
main function of the routine. See the example below. 


MEAD * main routine for processing a MEAD product. 
MEACC - read MEAN control cards. 

MEAINT - initialize MEAD variables and common 
blocks. 

MEAMTX - set up MEAD scaling matrix 
MEATRA - tremslate one line of input values into 
one line of output values. 


Exang>le 1 
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Han variables for constants. In U»o exani|>le below, constants 
such as Fortran unit numbers and the buffer nires are declared 
as variables. Su^ a convention facilitates pro9ram 
maintenance and revision. 


LOCAL VARIABLES 


REAL • 4 TIRE 

INTEGER • 4 BLANK/* •/, 
I hI*3F /Z3F/. 

I IKUNT/12/. 

i MAXLC/lOOCC/t 

i TRK7 /•7TRI'*/ 

LOGICAL * 4 ICPLG 

INTEGER * 2 LAROATI50C0I t 

LOGICAL « I lABUFClOOeit 


FUN 1/21 /t F20NT/22/t F3URT /23/t 
HEXFF /ZFF/t INBCKT /1008/» 
HAXCKK/3/f FAXIK/500/t 
NQ/*AC */t CUT10C200lt CUTLNT/11/* 


ROIL /Z7FFF/ 
ZERO/200/ 


RLiNK 

FIONT 

F2UNT 

F3UNT 

HEXFLG 

HEX3F 

HEXFF 

UeCNT 

imL 

MAXIN 

RAXLC 

OL'TUNT 

NC 

TRK7 

ZERO 


local VARIABLE DESCRIPTICNS 
ThE CCNS1ANT BLANK. 

CISK LNI1 WFLRE FIRST TAPE FILE IS TRANSFERRED. 

DISK LNI1 WHERE <FCCND TAPE FILE IS TRAKSFERREC. 

DISK LNtI WHERE THRIO TAPE FILE IS TRANSFERREC. 

EQUALS TIE IC FLAG FCR THE INPUT TAPE (DEPENDS CN 
FORMAT OF INPUT TAPE). 

CONSTANT EQUAL TO 3F HEXIOECtPAL. AN INPUT RECCRC IS 
AN ID RECORD IF THE FIRST BVIE ECUALS HEX 3F (7 TRACK FORMA 
CONSTANT EQUAL TO FF HEXOECIPAL. AN INPUT RECCRC IS 
AN IC RECORD IF THE FIRST BYTE ECLAIS HEX FF (9TRACK FORMAT 
NUMBER Of BYTES IN AN INPUT RECCRC. 

UNIT NUMfER OF INPUT TAPE. 

MAXIPLM NUMBER OF DATA CHANNELS THIS ROUTINE CAN HANDLE 

MAXIMUM NUMBER OF DATA VALUES IN CNE INPUT RECCRC. 

MAXIMUM number OF BYTES ALLCVEO IN CNE LINE CF LARSYS DATA. 
UNIT NUMBER FCR OUTPUT TAPE. 

CONSTANT EQUAL TC *NC*. 

CONSTANT EQUAL TO 'TTRK*. 

CONSTANT BYTE EQUAL TO CC HEAIDECIMAL. 


INPUT RECCRC IS 
HEX 3F (7 TRACK FORMAT 
INPUT RECCRC IS 
HEX FF (9TRACK FORMAT) 


Example 2 

In the above example, local variable descriptions have been provided 
only for the "constant" variables. See the example software system 
in Appendix C for descriptions of all local variables. 

4. Block commons must be named and they must have variables 
listed in the order: , g 

REAL * 4 
INTEGER * 4 
U>GICAL * 4 
INTEGER * 2 
UXJICAL * 1 





within each type of variable* the variables must be listed alphabetic-ally. 
Lar9e coronon blocks must be spaced for legibility. 


VAHUBLE NAMES fOU AGRONOMIC 10 AM) FOR SOUS 


QMMON 
QMMUN 
OMNON 
CQMM^ 
COMMON 
COMMON 
COMMON 
TOMMON 
^MMON 
COMMON 
COMMON 
MMON 
MMON 
UMMUN 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
.OMMON 
COMMON 
COMMON 


/IDNANC/ 
/ ONAMI/ 
/ DNAMC/ 

/IOnamI/ 
/ QNAMl/ 
/I INAMI/ 
/ ONANE/ 
/ OhAME/ 
/iONAME/ 
/IDNAME/ 
/IONAME/ 
/lOhANl/ 
/XDKANi/ 
/IONAME/ 
/IONAME/ 
/IONAMI/ 
/ ONANE/ 
/IONAME/ 
/{CNAME/ 
/IONAMI/ 
/IONAMI/ 
/ ONANE/ 
/IONAMI/ 
/iCNANl/ 
/IONAME/ 
/ CNAME/ 
/ ONAME/ 
/ ONAME/ 
/ CNAME/ 
/IONAME/ 


AIEE* 

CARO# 

C0MMO7I 

OBFR. 

OHSft 

0ENAI2)( 

OQF?l?r 

DQFTClli 

EPOAt 

EP09t 

LAlOt 

lF8t 

iSH: 

P^fSlA) 

RECAf 

ROUl, 

SENAU) 

TExfUl 

VI 

WABAIA) 

YEOA. 


A£IR« 

CAIN* 


OBGL 
DBTO t 
, OERA. 

• B8-”- 

EPlOt 

innUa) t 
loca(a) t 

SSfilii!; 

peSi' 

PLQAt 
t RATE 
REOA, 
RUSEt 
f SENO* 

TA TEt 

• TIOA, 

t WBTEt 
YELOt 


A/VI* 

CLCCt 

OAOA. 

OBNEf 
OIGRt 
Q0FAI2) • 
EPClt 


BAPRf BRLE 

CLTyIaI. C0B2 

OAFLt ceet 



EPC6, 


EXNAIA) 


FIVI* 


GMCS* 


HI SC* 

A) • 

iNsn 


LEPL* 

A) * 

LOOA* 

A) * 

MODA* 

A) » 

NMAT* 


OTST. . 


PLMCt 

REHOf 
SAGR, . 

specIa) • 

TANit 

TSNT, 

WEECf 

YELE* 


OBYL« 

DIIN, 

0GFSt2)f 

EP02* 

EP07, 

EXNU* 

FLLK2)t 

GRLCt 

H0RIC2I t 

JUOA 

LFIb).^ 

LOLAlilt 

MOFI lA) t 

NUDE* 

PECL* 

PHKO* 

PLNU* 

RENUt 

SCRA. 


OEEO 

DQFtC2) 

CCF6C 2) 

EPO'i 

EPOB 

FANAI Al 

FOCA 

HAWI 

ILLUC2) 

LF7 

L0L0I2) 

MOLA 

NUSA 

PE6R 

PHSE(A) 

PMOW 

SC?Vf A) 


STCOhOltSUCCI A) 
TCRl* TCR2 

UNCAt VARIIA) 


VARIABLE NAME* EXCLUSIVELY FOR SOILS ID 


COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 


/ICNAME/ 
/IONAME/ 
/IONAME/ 
/IONAME/ 
/IONAMI/ 
/ ONAME/ 
/IONAME/ 
/ ONAME/ 
/ OMAN I / 
/IONAMI/ 
/ ONAMI/ 
/IONAMI/ 
/IONAME/ 
/IONAMI/ 
/ ONAMI/ 
/IONAMI/ 
/ ONAME/ 


ACTI, 

BASA, 

E8SS: 

r.»; 

GRGR l2 1 1 

M|Uj 

MSNU* 

PASI* 

B: 

ii\r. 

WACO* 


ALUM* 

BUOE* 

CLAY* 

mi: 

hue!* 

LISH* 

MICL* 

OMOO* 

PHYS* 

SAPO* 

SOOI* 

SUBO* 

VALU* 

NAPH* 


ASH0I2) * 

BUFH* 

COCO* 

eSNU* 

EP13. 

HUE2* 

L0F161* 

MOTE* 

OKCA, 

PLIN* 

SHLl* 


WIOI* 

ZEIR* 


AVPH* 
CAEX* 
COIN* 
ELCO* 
ERFA * 

hsi* 

IRON* 

PAGN* 

MC2CI21* 

OROR* 

PLLI • 

SHRA* 


WISP 

ZEVl 


AVPC 

C/LC 

COPA 

ELNU 

EROS 

FSAN 

LlIN 

PANG 

PSAN 

PAPA 

PCTA 

SILI 


SOEL* SPGR* STAP 

SUCEIlOlvSUNAIA) » TEREI2I 
VeSA, VFSA* veSH 

WIER* YEAR 


INTERMEDIATE VAPIABLES LSEO IN CALCLLATICN CF IC VALUES 
PEAL * A VARI/BLES 


COMMON 

/IONAME/ 

ME TROW* 

XFRCO* 

XPLCO 



INTEGER 

; • A AARIABLES 





COMMON 

COMMON 

COMMON 

/{ONAME/ 

/IONAME/ 

/IONAME/ 

SURFST* 

BIOPLT* 

MOZ* 

SWING* 

FWR* 

ORC. 

TCPR, 

CRG. 
SI OE * 
TEX, 

HEAP 

SU6R 

INT8UH50J 

REAL • 

A PETFOM 

* XFRCO* 

XPLCO 






Although the common block above does not list the variables in the order 
REAL, INTEGER# LOGICAL# it is a good cxair^^le of spacing for legibility. 
(The variables are arranged by usage in this coirnnon block) . 
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Coiimon block variables must be describ(;d in a BLOCK DATA 
routific or in an initialization subroutine. The v.iriablc 
descriptions must be alphabetic. Sec Appendix D for an 
cxaii^le. 

5. Do not use Fortran entry points unless the use of them 
is clearly the best solution to an implementation prd^lem. 

€. Information and error messages should )>«■ informative to 
the user as well as the progrd'.ner. L.ieh message must 
include the name of the routine printiiiq the message. 



c 

•cLSE ERROR 


4200 

Ok • •FALSE* 

WRITE ITYFWTR. 942101 It ID(STRESSjn) 

WRITE IPRNTR# 942101 1* lOISTRLSSIin , . 

1 

94210 

FORMAT!* EOOAQ STRESS(*« 12* *) DOES NOT EQUAL** 



* N, Y* OR BLANK. *, 


2 

* INSTEAD IS SET TO I** A4* *)*• 

1 

3 

4X* *(XTKAI*1 


Exanf>le 4 


Zt may be numbered either sequentially (1 to n) or for the 

labeled Foreran statement nearest the message in the code. i 

In the exatr^le ed^ovc the message is numbered sequentially. 

In the exanple software system in Appendix C the messages 
are numbered for labeled Fortran statements. 

7. Labels for code statements must be assigned in ascending 
order within the body of each routine. For exaunplcs see 
i^pendix C. 

8. Labels for FORMAT statements must be assigned in ascending 
order within the body of each routine. The FORMAT label.s 

should be sufficieittly different from the code la'i' Is that | 

they stand out. For example, code labels in a routine could 
range from 100 to 900 and FORMAT labels from 9100 to 9900. 

The FORMAT statements may be intersperssd with the executable 
code or they may be just before the END statement. However, 
within one routine, they roust be either all inters|)ersed or 
all at the end. Tlie software system shown in Appendix C is 
an exait^le of FORMAT statements interspersed with executable 
code. 




j 
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9. Do not use unnecessary EQUIVAUBNCE statements. However, 
there are some data structures for which EQUIVAIENCE 
statements are necessary. For example, a LARSYS ID record 
contains real data values and integer data values. In order 
to correctly access both data types, the ID record must be 
declared as: 

REAL * 4 RID (2000) 

INTEGER * 4 ID (200) 

EQUIVALENCE (ID(1) , RID(l)) 

10. Use standard LARSYS and Reformatting routines whenever 
possible. For excui^le, often used LARSYS routines are 
CTLNRD and BCDVAL (for interpreting control cards), and 
often used Reformatting routines are IDRITE and EOT (for 
mounting LARSYS data tapes, writing ID records, and writing 
end-of-tape records) . 

11. Document all revisions to routines by adding your name and 
date to the ccTtiments. Include a version number if 
appropriate. If the revision is appropriate for only a 
special application, add a comment near the revision 
comment stating exactly what the special applications is. 


C 

C WRITTEN 07/19/79 BY CATHERINE KOZLOWSKI FOR FY70 

C SR&T CONTRACT 

C 

C REVISED 11/20/79 BY CATHERINE KOZLOWSKI FOR FY79 

C SR&T CCMJTRACT 

C 


Example 5 


12. Indent (horizontal) and space (vertical) the source eode 

to improve readability and/or logical flow of each routine. 
See the software system in Appendix C for examples. 



lnc« «•«««*««««*«*****«*«*#*«** ****««*«« 

AGRONCMIC RECORD SHEET NUMBER 2 CF THE GROUP OF 7 

RcAC (AGSHEt 9200* EN0>11C} 

AGUC2t PAGE2* RIOIHEIG). RIC(LEPU« IO(GRLE)t 
lOlYELEIt ICIBRLEI* XPLCQt XFRCG# 

METRQW, IDfPEGRI, RIOCOBGLl • RiDtORYUt 
PIO(CBBL)« RIDCCBST), RlD(DeFRK RIOIDBTC)« 
KIOIFRBD* BIOPLT* RID(LEAH)t IDCPLPC) 


F0RKATU3, II, 
FA.0,F3.0, 
FA.O, IX, 12 


F3.0, IX, FA.l, 


2F6.0, 


F5*0, 


5F5.0, 


PTRA s PTRA - 1 
KEADIAGSHE, 9201, ENOsllO) 
BLKIOfFEIG), RLKIOCLEPL), 
BLKID(eRLE), BLKID(PLCO), 
BLKlDfCBGLI, BLKID(DBYL), 
BLKICICBFR), BLKIO(OBTO), 
BLKIO(PLMO) 


BLKIDfGRLE), BLKIO(YELE), 
eLKtO(FRCL), BLKIO(PEGR), 
BLKIO(DBBL), BLKIO(DBST), 
BLKID(FRBl), BLKID(LEAR), 


BLKlOlPLMO} 

F0RMAT(T5, A3, IX, A2 , TIA, 3/>2, IX, 2A3, T32, 
IX, 511K,AA), 2(2X,AA), T76, A2 ♦ IX, A2) 


Example 6 


If possible, use the following convention for FILEDEFing 
and assigning tape units: 

FILEDEF 11 TAPI 
FILEDEF 12 TAP2 
FILEDEF 13 TAP3 
FILEDEF 14 TAP4 
FILEDEF 10 TAPS 


where Fortran unit 11 is the output tape and units 12-14, 
10 are input tapes. 
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15. Several suggestions about labels and CONTINUE statements: 

a. It is easier to revise routine if each DO loop has 
its own CONTINUE statement. 


DO 120 K > 1, 20 
DO 100 J = 1, 3 

ARRAY (J,K) = J + K 
100 CONTINUE 
126 CONTINUE 


Example 7 

b. It is easier to revise a routine if all of its non-> 
FORMAT labels are on CONTINUE statements. 


16 . Debugging convention 
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B. Layout of Individual Routines 

C routine name 
C 

Q******************** ************************* ***************** 

C routine name one- line description 
C 

C WRITTEN date BY name FOR CONTRACT name or number 
C REVISED data BY name 
C 

StBROUTINE name 
C 

IMPLICIT INTEGER *4 (A - Z) 

C 

C detailed description 
C 

C special features and/or limitations 
C 

C input 
C 

C output 
C 

C subroutines used (include one-line description of 
each subroutine) 

C 

COMMON /name/ declarations 


COMMON / name/ declarations 
C 

0 **** ********************************************** ************ 

c 

C LOCAL VARIABLES 

C 

local variables declarations by type, then alphabetic 
(include parameters as necessary) 

C 

C local variable descriptions including parameters, listed 
alphabetically 
C 

body of routine 
END 


Example 8 


All routines should follow the general format outlined above « 
See Appendix C for a complete system following this layout. 





ooonoooooo ooooo onooonor>oor>or>r>oooo| I 
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1. The first several lines of the source code should identify 
the routine. .< 


SP.CSC» 


SPCSCN REFCRMfTS CNE SPECSCAN TAPE TO CNE lARSVS RUN. 


HPITTEN 02/l«/7< BY CATPERINE KOUOESKI 
REVISED CA/C4/7< BY CATHERINE KCZtCVSKI 
REVISED C7/0</7S BY CATHERINE KOZLOVSKI 


FCR SRtT FY79 CCNTfACT 


Example 9 


2. After the IMPLICIT INTEGER * 4 statement, there should be 
a detadled. description of the routine . 


SPCSCN 

, THIS PROGRAM AND ITS SUERCUTINES REFCRMAT A 7-TRACK (MCCE 31 
OK 9-TRACK 8CC BPI SPECSCAN TAPE TO A 9>TRACK 1600 BPI LARSYS 
OATA RUN. THE CRIGINAL SPECSCAN TAPE FCR NHICH THIS SCFTKARE WAS 
WRITTEN NAS RECilVEC FRCP RC6ERT A. GGOCOIN'G* TECHNICCLCR GRAPHIC 
SERVICES INC. {TEXAS OPERATICNS. LYNOCN B. JCHNSCN SPACE CENTEK» 
P.O. BOX 56863* HOUSTON, TEXAS 77C56I. 

THE INPUT TAPE HAS CNE OR PCRE FILES, EACH FILE CCRRESPCNCINC TO 1 
CHANNEL OF CATA ON THE LARSYS TAPE. ALL RECCRCS CN THE INPUT TAPE 
ARE 1C08 BYTES lONG. THE FIRST INPLT RECORD OF EACH FILE PAV 
BE AN IC RECCRO -- IF IT IS, THE FIRST BYTE ECUALS HEXICECIMAL 
T3F* OR 'FF* ANE THE RECORD, I S SKIPPED DURING PROCESSING CF CATA. 
IT IS ASSUMED EACH FILE HAS THE SAME NUHEER OF RECORDS AND, 

IF ONE FILE HAS AN ID, THEY ALL HAVE ID RECCROS. 

IT IS ALSO ASSUMED THAT THE ONLY SIZE OF CNE SPECSCAN INPUT 
RF.CORC IS 1CC8 BYTES. HOWEVER, ONE LINE OF SPECSCAN INPUT MAY 
BE SEVERAL INPUT RECORDS LONG. 

CONTINLE 

THE PROGRAM REQIIRES ONE TEMPORARY CISK ANC TWO TAPE CRIVES. 

THE EXEC CACLLATES THE AMOLNT CF TEMP SPACE NEEDED FOR THE DATA 
TC BE TRANSFERRED. THE INPLT TAPE DRIVE IS ASSIGNED UNIT NUMBER l: 
IT MAY BE A 7 OR 9 TRACK DRIVE. THE OUTPUT TAPE DRIVE IS ASSIGNED 
UNIT NUMBER 21 IT IS A 9 TRACK DRIVE. 

CCNTIME 

THE PROGRAM FIRST TRANSFERS THE INPLT TAPE FILES TC DISK. THEN 
IT SETS UP THE lARSYS ID RECORD. INPUT OATA IS REFCRMATTED LINE 
BY LINE. EACH INPUT LINE GCES THRCLGH THE FCLLCWINC PROCESS INGO 
1) INPUT FUF CHANNEL N IS READ FROM CISK N INTO A TEMP BLFFER. 
21 THE CATA IN THE TEMP BUFFER IS REFORMATTED FRCM 2 BYTES 

PER VALUE TO I BYTE PER VALUE. THE CATA IS INVERTED (BLACK 
TO white: at this TIME. IF NECESSARY, THE CATA IS ALSO 
FLIPPED LEFT TO RIGHT. 

31 THE CATA IS MOVED TC A LARSYS-FORMAT CATA LINE 
AND WRITEN TO TAPE. 


TWO TAPE CRIVES. 

NEEDED FOR THE DATA 
ASSIGNED UNIT NUMBER 12 
TAPE DRIVE IS ASSIGNED 


Examole 10 


In the cibove example, special features and limitations of the 
routine have been noted. Special features are 1) the input can 
be on either a 7-track or 9-track tape, and 2) the data can be 
flipped left to right. Limitations are 1) if one input file has 
an id record, all input files must have ids, and 2) the routine 
requires two tape drives and one tenporary disk. 
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3. Input requirements must be specified. 


C 

c 

c 

c 

c 

c 


THE INPUT IS AN AFRAV OF DATA VAltES IN SFECSCAK TAPE 
FORMAT. 

EACH CATA VALUE IS ASSUMED TO BE CNE 3>BI1 FIELD IN 
AN 8-6IT BYTE AND ONE 6-BIT FIELD IN AN 8-BIT BYTEt 
THESE TWO FIELDS REPRESENTING CNE DATA VALUE RANGING 
FROM 0 TO 511. 


Exan^le 11 


4. Output from a routine must be described. 


C 

C 

c 

c 


THE OITPUT i; AN ARRAY CF DATA VALUES WITH EACH 2 BYTES 
REPRESENTING ONE 8-BIT CATA VALUE ITI-E FIRST BYTE IS SET TO ZERC 
AND THE SECOND BYTE CONTAINS THE CATA}. CUTPLT VALUES RANGE 
BETWEEN 0 AKC 255. 


Example 12 


5. The source listing must include all non-system stibroutines 
called. 


c 

THE NON-SYSTEM 

SUBROUTINES USED AKEC 

c 

EOT 

VRITES END-CF-TAPE RECCRC N CUTPUT TAPE. 

c 

GTDATE 

RETURNS TOOAYS'S DATE IN CNARACTER FORMAT. 

c 

IDRITE 

►OUNTS OUTPUT TAPE AND WRITES CUTPUT RUN 10 RECORD. 

C 

MOUNT 

NOUNTS INPUT TAPE. 

c 

MCVBYT 

MOVES BYTES FROM INPUT BUFFER TC CUTPUT BUFFER. 

c 

SPCDAT 

TRANSLATES CATA FRCM A 9-BIT FORMAT TO AN 8-BIT 

c 


FORMAT. 

c 

SPCSAP 

CALCULATES THE NUMBER OF SAMPLES PER CHANNEL IN THE 

c 


OUTPUT. 

c 

TAPOP (ENTRY POINTS TQPEF. TCPFF, 1 

C 

c 


TOPRO, TOPV.R) PERFCRM.S TAPE I /C FUNCTICNS. I 

c 

C 

c 


Example 13 



oo«->r>r>r»rvoor>ooo( 
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6. All local variables must be declared (as necessary) and 
described. 


C**** '>*■***** ***«<■** ******** «*«« 4 4* a* 4 sees *««<#«**«•* 

C 

INTEGER « 4 BLFCNT /5C4/« NAXOAT /500/t ACAOAT /4/» RLIFIT /5/t 
^ 1 ZERO 70/ 

^ INTEGER * 2 BIFFERISOA) 

^ LOGICAL ♦ 4 IlFLG 

^ LOGICAL « 1 UeOFdOCEl 

^ equivalence (EUFFERt INBLFI 

C****44************4*******************4***4****4«**4*****4«**4*4*4***** 

local variable DESCRtPTIGNS 

ACJUST TEMPORARY VARIABLE LSEf) TO AC JUST SAFPLE CCUNT TC BE EVENLY 
OIVISIBLf BY 4. 

BUFCNT NUMBER OF ELEMENTS IN INPUT BUFFER. 

BUFFER INPUT BUFFER IN INTEGER * 2 FORMAT. 

CISP CISPLACEPENT INTC INPUT BUFFER. 

ICFLG FLAG INDICATIANG WHETHER FIRST INPUT RECORD CF THE FILE IS 
AN ID RECORD. IF ICFLG IS SET# FIRST RECORD IS AN ID. 

INEUF INPUT BUFFER IN LDGICAL * I FORMAT. 

LSTVAL LAST CATA VALUE IN INPUT BUFFER. 

MAXDAT MAXIMUM NUMBER OF DATA VALUES POSSIBLE IN CNE INPUT RECORD. 

CONTINUE 

C NCNCAT NUMBER OF NON-DATA VALLES IN INPUT BUFFER. 

C NREAD NUMBER OF CONSECUTIVE RECORDS READ WHEN SEARCH FCR ZERO DATA 

C VALUES. 

C NSAMP NUMBER OF SAMPLES PER CHANNEL THAT WILL BE IN LARSYS CUTPUT. 

C PPEVAL PFEVICUS DATA VALUE IN INPUT BUFFER. USED TC SEARCH 

C BACKWARDS IN INPUT RECORD. 

C RLIMIT UPPER LI FIT CN THE NUMBER OF CONSECUTIVE READS TC PERFORM 

C BEFORE TIRMINATING SEARCH FOR ZERC DATA VALUES. 

C TOTAL RUNNING TOTAL USED TC CALCULATE NUMBER CF SAMPLES. 

C UNIT DISK UNIT FROM WHICH TO READ INPUT RECORDS. 

C ZERO THE CCNSTANT ZERO. 

C 

C444444444«*****4**4****4*4***4********44*«44*444*4»4*4**4****4*4******* 
C444****4«****4*444*444******* 4******4*44**44*4*444*444*44**444*4*4*4*4* 

C 


Exan^le 14 
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C. Comments Within Ihe Body of a Routine 

1. Highlight comments that describe large sections of code. 

See Appendix C for examples. 

2. Comments by themselves should describe thr flow of the 
routine in sufficient detail so a reader can understand the 
routine without looking at the code. 

3. Inobvious programming "tricks” must be explained in detail 
including the reason for the trick. 

4. Specific suggestions: 

a. Comment a control card confuted GO TO so that it is 
apparent which label corresponds to which key word. 


r 

IMPLICIT IKTECER ♦ 4 

(A - ZI 


r 

INTECifR ♦ 4 

KEYLST(7) 




L'ATA KEYLSl 

/ •*INP«, 

•RcFU*. 

•INPL't *SCRA*, *CUTP*, 

r 

L 

•ENDS 

•-COM*/ 


V 

c 

r 

DATA KEYSZ 

n ! 



L 

c 

r 

CALL CTLWRC(C/RD, COL 

f KEYLST 

t KEYSZ* CODE* REAUIN* ERRCR) 

C 

♦ IMP 

REFO 

INPL 

SCRA CUTP EhD -CCF 


GOTO (lOOOt 

1000* 

300Ct 

4000* 5000* 6000* 7000)* CODE ! 

looa 

CONTINLE 




2000 

CONTINUE 




30C0 

CONTINLE 




4000 

CONTINUE 




5000 

CONTINUE 




6000 

CONTINLE 




7000 

CONTINUE 





STOP 





END 





Example 15 




non non 
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Comment logical program structures with statements 
such as: 


WHILE NOT END-OF-FILE PROCESS DATA 


REPEAT LINE PROCESSING UNTIL END-OF-FILE 


C 

C IF GOOD DATA THEN PROCESS IT 
C ELSE PRINT ERROR MESSAGE 

C 


APPENDIX A 


07-20-79 


A4 


* * • CFMRP Control Cards * • • 


Keyworu 

R 

E 

2 

Option 

Parameter 

Function 

Default 

INPUT 


TAPE 

(xxxx) 

Input tape number. 

(none) 

OUTPUT 

* 

TAPE 

(yyyy) 

Output tape number. 

(none) 



FIRST 


Indicates that output 
run sequencing starts at 
beginning of tape.^ 

(none) 



LAST 


Indicates that output run 
sequencing starts after the 
last run sequence already 
on tape.^ 

(none) 



RUSE 

(aaaa) 

Start output run sequencing 
at output run sequence 
number aaaa.^ 

(none) 

INRUS 


FRUSE 

(bbbb) 

Indicates first input run 
to be CFMRPed.^ 

bbbb- 'FIRS’ 



UlUSE 

(cccc) 

Indicates last input run 
sequence to be CFMRPed.^ 

cccc-999999999 

END 

* 

(none) 


Indicates end of a processing 
group. 


SEND 


(none) 


Indicates end of all 



processing. 


NOTES; 

1. Must have a FIRST, LAST, or RUSE option. 

2. The defaults are set up so the entire input tape would be CFMRPcd. Therefore, 
a deck setup of 

INPUT TAPE (9900) 

OUTPUT TAPE (9901) , FIRST 

END 

SEND 

would CFMRP all the run sequences on tape 9900 anu </Ut the output on tape 9901 
starting at run sequence 1. 


45 


07 - 20 - 7 *> 


* • * CFMRP Control Cardb * • * (coni) 


A deck setup of 

INPUT TAPE(9902) 

OUTPUT TAPE (9903), RUSE (S) 

XNRUSE FRUSE(7), LRUSE(IS) 

INRUSE FRUSE(20), LRUSE(2S) 

END 

SEND 

would CFMRP run sequences 7 to 15 inclusive and 20-25 inclusive on tape 9902 
and put the output on tape 9903 as run sequences 5 to 19 inclusive. 
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APPENDIX B 


L47S 


, LARS Program Abstract 
47 

module IDEWTinCATION 

Module Name: IDEOT Function Name:^ 


Purpose: Edit ID record and modify data of MSS tape 


System/Language t 360 /fortran 



Author: Laura Wallace 

Date: 7 / 28/78 


Latest Revisor: Laura Wallace . _ 

Date: io/05/78 



MODULE ABSTRACT 


Through eoatrel cards, the ZD record of a LARSYS formatted MSS data tape 
(either 800 or 1600 BPZ) la edited and output to another LARSYS formatted 
tape (1600 BFl). More than one input tape may be concatenated onto one output 
tape if the total line count conprlaea leas than 9SZ of the output tape. There 
la an option to remove ancillary data In the data records (see MOOTAP subroutine 
abstract). There Is also an option to ROTATE the data 180 degrees, north-south 
(see ROTDAT subroutine abstract). 


UNIVERSITY 

LaboratT/ry for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 


1. MODULE ABSTR ACT! 
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Input is expected from the card reader (device 5). 

Mote t The FILEDEF in the EXEC allows input to be on disk. 

B. EXEC Cards Needed 

GETDISK LARSYS ADR 19A MODE A 
LOAD IDEDT (CLEAR NOMAP) 

USE CTLVRD BCDVAL CPFUNC 
FILEDEF 5 DSK-Pl & 1 & 2 
START IDEDT 
FILEDEF * CLEAR 
& EXIT 

C. Running IDEDT from the Terminal 

IPL REFORM 

IDEDT filename fUetype (RUN IDEDT FROM ABOVE EXEC) 

D. Running IDEDT from Batch 

Batch Cards ; 

BATCH MACHINE BATONITE 
BATCH ID userid username 
BATCH OUTPUT printloc punchloc 
EXEC$$ 

GETDISK REFORM 19D D 

GETDISK userid 191 P RR DETACH PASS diskpassword 
GLOBAL T REFRMLIB SYSLIB 
EXEC IDEDT filename filetype 
$$ 


2. INTERNAL DESCRIPTION 


All control cards are read. ’ OPTION MMS * puts spectral band maximums and 
minimums into a temporary ID array (TEMPID) . EDIT puts new information into 
the same TEMPID . When the * END * is read all processing begins. The first 
input tape is mounted and the ID is read. The number of lines in each of the 
input tapes is counted and printed. The first N tapes are concatenated onto 
the output tape. N is the number of tapes whose sum of total line count comprise 
less than 95Z of the output tape The ID record is then read from the input 
tape and edited using the TEMPID . IDRITE gets the output tape and writes out 
the ID. If the data is to be rotated 180 degrees then the tape is forward 
spaced the number of lines to be written. Each line of data is then read (using 
TOPRB if rotate option selected) from the input tape, modified if the MMS option 
is given, rotated if requested, and written out to the output tape. The 
LARS17 forms are printed. The tapes are rewound and detached. 
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SUBROUTINES CALLED 
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BCOFIL 

MOUNT 

CPFUNC 

RCOUNT 

CTLPRM 

RINGIN 

CTLHRD 

ROTOAT 

EOT 

TOPBF 

IVAL 

TOPFF 

GTMTE 

TOPFS 

lORITE 

TOPRB 

IVAL 

TOPRI> 

LARS17 

TOPRU 

MODTAP 

TOPWR 


4. INPUT DESCRIPTION 


E 

9 -. 

RET 

WORD(COL.l) 

CONTROL 

PARAMETER 

FUNCTION 

DEFAULT 

* 

INPUT 

TAPE(XXX) 

MAX OF 5 INPUT TAPES 

00000 



FILE(FF) 

LOCATION OF FILES FOR INPUT 

1 



RUN(XXXXXEXX) 

TO BE IMPLEMENTED AT A LATER DATE 

(NONE) 



DENSITY (XXXX) 

DENSITY OF INPUT TAPE 

800 BPI 


OUTPUT 

TAPE(XXX) 

OUTPUT TAPE (1600 BPI) 

00000 



FILE(FF) 

FILE LOCATION FOR OUTPUT 

00000 



RUN(XXXXXXXX) 

RUN NUMBER FOR OUTPUT TAPE 

(NONE) 


OPTIONS 

MMS 

MODIFY DATA, ADD SPECTRAL BAND 

(NONE) 




INFOR FOR MODULAR MULTISPECTRAL 





SCANNER (MSS) 




FORH(X) 

NUMBER OF IARS17 FORMS TO PRINT 

8 



DEBUG 

PRINT ADDITIONAL INFORMATION AND 

OFF 




VARIABLE VALUES 




ROTATE 

ROTATE DATA 180 DEGREES. ONLY VSBBB 

DO NOT ROTATE 




VALID FOR 1 INPUT TAPE. IDfl6) 





AUTOMATICALLY CHANGES BY 180 





DEGREES. 




NOCOUNT 

DO NOT COUNT THE RECORDS ON THE 

COUNTS THEM 




TAPE 



LARS17 

MISSION (XXX) 

MISSION NUMBER 

0 



SITE (XXX) 

SITE NUMBER 

0 



LINE(X) 

LINE NUMBER OF LINE-RUN- 

0 



RUN(X) 

RUN NUMBER OF LINE-RUN- 

0 

* 

EDIT 

X 

ELEMENT OF ID TO BE EDITED 

(NONE) 



NNNNN 

NE17 VALUE (CAN BE INTEGER, REAL, 

(NONE) 




OR ALPHANUMERIC) 


* 

END 


BEGIN PROCESSING 


NOTE: The run 

number for the output tape is the same run number as in 

the input 


tape id 

record unless the 

'RUN' parameter is specified on the '( 

output ' 


keyword or 'edit* for ID(3) is used. Since the keywords are operated on 
in order, the keyword which is later in the list of keywords input will 
take precedence over previous changes. 


5. OUTPUT DESCRIPTIOM 
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Control cards vlll be printed, followed by the number of lines of data on 
the input tapes (unless NOCOUNI option specified) and which tapes were actually 
read and the new value for the number of data lines, as pul in ID(20) . 

The previous values and new values for all edited elements of the ID record are 
printed. IDRITE prints the percentage the input will take and any information 
pertaining to the run number, followed by the IDPRIHT of the tape. The IARS17 
f oras are printed (caused by LARS17 subroutine) . 

ERROR MESSAGES: (All errors cause the program to halt. The tapes are rewound 

and detached.) 

ERROR OH CARP — — (the card is echoed) 

JOB HALTED No output written on tape 

INPUT TAPE NOT WRITE PROTECTED . HALTED . 

** EOT DETECTED ON INPUT TAPE BEFORE ALL LINES READ 
** WITH NO 2ND INPUT TAPE GIVEN 
** NUMBER OF LINES ACTUALLY WRITTEN - XXXXX 
**N0 MORE ROOM ON OUTPUT TAPE. LINES WRITTEN = XXX 

***RUN NUMBER FOR OUTPUT TAPE IS ILL DEFINED — JOB TERMINATED 
ERROR DURING TOPRB. PROGRAM TERMINATED. OUTPUT TAPE PARTIAtTy 
WRITTEN THRU XXX RECORDS . 

ILLEGAL CANNOT HAVE MORE THAN ONE INPUT TAPE WITH ROTATE 

OPTION JOB TERMINATED. NO OUTPUT . 


6. SUraOUTINE AND ENTRY POINTS CALLED 

RCOUNT 

ROTDAT 

CTLWRD 

CTLWRD 

CTLPRM 

BCDVAL 

IVAL 

FVAL 

BCDFIL 

TAPOP 

TOPHR 

TOPRD 

TOPRU 

TOPBF 

EOT 

CPFUNC 

IDRITE 

LARS17 

MODTAP 

REVERS 

7. SAMPLE CONTROL CARD DECK 

INPUT TAPE (26,27), FILE(1,1), DENSITY (800) 
OPTIONS MHS,F0RM(4) 

LARS17 MISSION(365) ,SITE(194) ,LINE(2) ,RUN(1) 

EDIT 3,77007500 

EDIT 7, LINE 2 RUNl HAND 

EDIT 21, 44.58 

END 



#«• en €>!•»• 


FIIEO SPCSCN EXEC 
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PURCUE / LARS 3931 


kPlTTEN 07/C2/79 PY CATPEMtAE ROZLCMSKI FCR SRCT FV79 CCNTRACT 
BEGTYPE 

POM PANY FILES hlLL BE TRANSFERRED TC DISK* HOk FANY RECCPOS 
PER FILE 

END 

CREAC VARS £KF ILE CNREC 

CALCULATE PCk PANY CYLINDERS Of 3350 TEPP SPACE IS AEECED 

CTOT s CAFILE « CNREC « ICOf 
CNCYL « CTOT / A60C00 
CNCYL * CKCYL ♦ I 

CTYPE CNCYL CYL OF 3350 TEMP SPACE kILL BE REOLESTO 

OEF T3350 25C CNCYL 

CIF CRETCCCE NE 91 CGOTO -GOTDSK 

CTYPE TEPP CISK WITH CNCYL OF SPACE ACT AVAIL. 

-GOTDSK CiPcRETCCCE ED 92 C60TC -fORMATD 
FCRMAT 250 E 

-FORMATD ACCESS 250 E (ERASE 

GLOBAL TXTLie REFRMLIB CMSLII FORTRAN 
LOAD SPCSCN TAPCP ( NOMAP ClEAR 
FILEDEF 5 READER 

6 PRINTER (RECFM FA 

7 PUNCF 
' TERP (PERM 
TERP (PERM 


15 

16 


FILEDEF 
FILEDEF 
FILEDEF 
FILEDEF 
OUTPUT TAPE 

FILEDEF II TAPI (PERM 
INPUT TAPE 

FILEDEF 12 TAP2 (PERM 7TRACK DEN 800 TRTCH C 
TEMPORARY DISK FILES 

FILEDEF 21 CISK SPEC FILEl B] (RECFM F 

FILEDEF 22 CISK SPEC FILE/ B1 (RECFM F 

FILEDEF 22 DISK SPEC FILE3 B1 (RECFM F 

START SPCSCN 
CP CET 181 
CEXIT 


LRECL ICOE BLKSIZE 1008 
LRECL ICOE BLKSIZE 1008 
LRECL 1C08 BLKSIZE 1008 
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PUROUE / LARS 3031 


SPCSCN REFCRMATS ONE SPECSCAN TAPE TO 

WHITTEN 02/1A/79 BY CATHERINE KOZLOtiSKI 
REVISED CA/OA/79 BY CATHERINE K02L0VSKI 
REVISED (7/02/7S BY CATHERINE KOZLOkSKI 


LARSVS 


FY79 CONTRACT 


IMPLICIT INTECER « A (A>2I 
SFCSCN 

THIS PROGRAM AND ITS SUBROUTINES REFCRMA1 A 7-TRACK IMCCE 3) 

OR 9-TRACK 8CC BPI SPECSCAN TAPE TO A 9-TRACK 1600 BPl LARSYS 
DATA RUN. TEE CRICINAL SPECSCAN TAFE FCR WHICH THIS SOFTWARE MAS 
WRITTEN WAS RECfiVEC FROM ROBERT A. GOCCDINC* TECHNICCLCR GRAPHIC 
SERVICES INC. (TEXAS OPERATIONSt LYNDON B. JOHNSCK SPACE CENTER* 

P.O. BOX 588(3* HOUSTON. TEXAS 770S8I. , 

THE INPUT TAPE FAS CNE CR PORE FILES, EACH FILE CCRRESPCNC INC TO 1 
CHANNEL OF DATA ON THE LARSYS TAPE. ALL KECCROS ON THE INPUT TAPE 
ARE 1008 BYTES LONG. THE FIRST INPUT RECORD OF EACH FILE PAY 


BE AN ID RECORD — 


the FIRST 


ECUALS HEXICECIMAL 


TWC TAPE DRIVES. 

NEEDED FOR THE DATA 
ASSIGNED UNIT NUNBER 12 
TAPE DRIVE IS ASSIGNED 


•3F* OR *FF« ANC THE RECORC IS SKIPPED CURING PROCESSING OF DATA. 

IT IS ASSUMEC EACH FILE HAS THE SAME NUMBER OF RECORDS AND* 

IF ONE FILE FAS AN 10* THEY ALL HAVE ID RECCROS. 

IT IS ALSO ASSURED THAT THE ONLY SIZE OF ONE SPECSCAN INPUT 
RECORD IS lOCe BYTES. HOWEVER, ONE LINE OF SPECSCAN INPUT MAY 
BE SEVERAL INPUT RECORDS LCNG. 

CONTINUE 

THE PROGRAM REQUIRES ONE TEMPORARY DISK ANC TWC TAPE DRIVES. 

THE EXcC CACULATES THE AMOUNT OF TEMP SPACE NEEDED FOR THE DATA 
TC BE TRANSFERRED. THE INPUT TAPE DRIVE IS ASSIGNED UNIT NUMBER i: 
IT MAY BE A 1 OR 9 TRACK DRIVE. THE CUTFUT TAPE DRIVE IS ASSIGNED 
UNIT NUMBER 11 IT IS A 9 TRACK DRIVE. 

CONTINUE 

THE PROGRAM FIRST TRANSFERS THE INPUT TAPE FILES TC DISK. THEN 
IT SETS UP TFE LARSYS ID RECORD. INPUT DATA IS REFCRPATTED LINE 
BY LINE. EACH INPUT LINE GOES THROUGH THE FOLLOWING PROCESS INGO 
1) INPUT FOR CHANNEL N IS READ FROM DISK N INTO A TEMP BUFFER. 
21 THE DATA IN THE TEMP BUFFER IS REFORMATTED FRCP 2 BYTES 

PER VALUE TO 1 BYTE PER VALUE. THE DATA IS INVERTED IBLACK 
TC WHITE! AT THIS TIME. IF NECESSARY* THE DATA IS ALSO 
FLIPPED LEFT TO RIGHT. 

3) THE DATA IS MOVED TC A LARSYS-FORMAT DATA LINE 
AND WRITEN TO TAPE. 

CONTINUE 

THE NON-SYSTEM SUBROUTINES USED AREC 

EOT WRITES END-CF-TAPE RECORC N CUTPUT TAPE. 

GTCATE RETURNS TODAYS'S DATE IN CHARACTER FCRPAT. 

ICRITE MOUNTS OUTPUT TAPE AND WRITES OUTPUT RUN ID RECORD. 
MOUNT FOUNTS INPUT TAPE. 

MOVBVl MOVES BYTES FROM INPUT BUFFER TO CUTPUT BUFFER. 
SPCDAT TRANSLATES DATA FROM A 9-BIT FORMAT TC AN 8-BIT 
FORMAT. 

SPCSAM CALCULATES 1HE NUMBER OF SAMPLES PER CHANNEL IN THE 
CUTPUT. 

TAPOP (ENTRY POINTS TOPEF, TCPFF , 

TOPRO* 10PWR) PERFCRMS TAPE I/C FUNCTICNS. 


THE NON-SYSTEM 
EOT 

GTCATE 

ICRITE 

MOUNT 

MOVBYI 

SPCDAT 

SPCSAM 

TAPOP ( 


l«««4«*)»*«****A«*«4*«***4>****»4>*4i4>A***4«*«4*«**4*«*««**«4*«4«4«4*4**4«* 


COMMON /SPCCOM/ 

INTECER « 2 VARIABLES 

1 INOAI* CUTOAT, 

LCGKAL * 1 VARIABLES 

2 FLIP 
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INTEGER » 2 IPOATIIOOCOI* 0UT0A1I lOOOCl 
LOGICAL * 1 FLIP 


LOCAL VARIABLES 


INTEGER 


LOGICAL 

INTEGER 

LOGICAL 


BLANK/* •/• 
HEX3F /23F/, 
IhUNT/12/* 
M/XLC/IOCCO/* 
TPK7 /*7TRK*/ 

ICFLG 

L/R0AT(5000I • 
lABUFllOOei* 


FlUNT/21/t 
HEXFF /2FF/« 
HAXCFK/3/, 
NO/*AC •/. 


ROLL /Z7FFF/ 
2ERO/ZOO/ 


F2UNT/22/t F3UM /23/, 
INBCNT /lOOe/t 
PAXIK/500/t 

CUTI0C200U CUTLNT/ll/t 


C*«*)l>«<t>«**«****««*>l>>»**>**««****«******«*4«**4*«*«4«*««*«««4***««*4***4>*<»* 

c 

C LOCAL VARIABLE OESCRIPTICNS 


CFILE 

DISP 

FLIP 


FlUNT 

F2UNT 

F3UNT 


QUESTIONS ASKEO 


TERPINAL. 


ANSM USEO TO STORE ANSWER TO QUESTIONS ASKEO AT THE TERPINAL. 

BL/NK THE CCNS1ANT BLANK. 

CHKIC FIRST BYTE OF FIRST RECORD OF A FILE. 

CCUNT BYTE COURT USED IN CALLS TO 10PRC ARC TCPWR. 

CURCHN CHANNEL CURRENTLY IK PROCESS. 

CURLIN LINE NUMBER OF LARSYS LINE CURRENTLY IN PROCESS. 

CFILE UNIT NUMBER OF DISK FILE CURRENTLY BEING PROCESSED. 

DISP DISPLACEPENT INTO LARSYS LINE CURRENTLY BEING ASSEPBLED. 

FLIP FLAG INCICATING WHETHER INPUT DATA NEEDS TO BE FLIPPEE 

LEFT TO RIGHT ON OUTPUT. IF FLIP IS SET, THE DATA SHCULO 
BE FLIPPED. 

CCNTINUE 

FlUNT DISK UNIT WHERE FIRST TAPE FILE IS TRANSFERRED. 

F2UNT DISK LNI1 WHERE SECCNO TAPE FILE IS TRANSFERREC. 

F3UNT DISK LNI1 WHERE THRID TAPE FILE IS TRANSFERRED. 

HEXFLG EQUALS THE ID FLAG FOR THE INPUT TAPE lOEPENDS ON 
FORMAT or INPUT TAPE). 

HEX3F CONSTANT EQUAL TO 3F HEXIDECIPAL. AN INPUT RECORD IS 

AN ID RECORD IF THE FIRST BYTE EQUALS HEX 3F 17 TRACK FORMAT) 
HEXFF CONSTANT EQUAL TO FF HEXDECIHAL. AN INPUT RECORD IS 

AN ID RECORD IF THE FIRST BYTE EQUALS HEX FF (9TRACK FORMAT) 
ICFLG FLAG INICCATING WHETHER FIRST INPUT RECORD OF A FILE IS AN ID 


HEXFLG EQUALS 
FORMAT 
HEX3F CCNSTA) 
AN ID I 
HEXFF CONSTA) 


ICFLG 


IF lOFLG IS SET. FIRST RECORD IS AN 
INBCNT NUMBER OF BYTES IN AN INPUT RECORD. 


FORMAT) 
IS AN ID 


NCAL 


THE VALLES 


THERE 


SPECSCAN 


TOTAL 


INEUF BUFFER FOR INPUT RECORDS. 

INCAT INTERPECIATE LINE BUFFER WITH THE VALUES IN SPECSCAN 
ITHAT ISt INPUT) FORMAT. 

CONTINUE 

INFIL FILE NUMBER OF FIRST INPUT FILE. THERE IS A TOTAL CF 3 
INPUT FllES. 

INTAP INPUT TAFE NUMBER. 

INUNT UNIT NUMBER OF INPUT TAPE. 

LAROAT BUFFER FCR A PROCESSED LINE CF LARSYS DATA. 

LLNT BYTE COUNT FOR ONE LINE OF LARSYS DATA. 

MAXCHN MAXIMUM NUMBER OF DATA CHANNELS THIS ROUTINE CAN HANDLE 
MAXIN MAXIMUM NUMBER OF DATA VALLES IN ONE INPUT RECORD. 

MAXLC MAXIMUM NUMBER OF BYTES ALLOWED IN CKE LINE CF LARSYS DATA. 
MCCUNT NUMBER Of BYTES TO TRANSFER FROM INPLT BEFFER TO 
INTERMEDIATE LARSYS LINE BUFFER. 

NCAL NUMBER OF CALIRRATICN BYTES FER CHANNEL IN INTERPECIATE 
LARSYS LINE BUFFER. 
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NCHAN 

NCAT 

NLINE 

NC 


NUMBER Of CHANNELS CF DATA. 

NUMBER Of DATA BYTES PER CHARNEL If 
L INE eUFFER. 

TOTAL NUPBER OF LARSVS DATA LINES. 
CONSTANT EQUAL TC *KC*. 


IRTERPEOIATE LARSYS 


NSAMP NUMBER Of SAMPLES PEfi CHANNEL PEP LINE f INCLUDING CALI6RATICN 
ONFOUT NUMBER Of INPUT RECORDS THAT CORPESPCNO TC ONE CHANNELt ONE 
LINE CF CUTPUT. 

OtTCAT INTERPECIATE LARSYS DATA LINE BUFFER. EACH DATA VALUE TAKES 
2 BYTES eUT ACTUAL GATA IS IN SECOND BYTE. 

OUTIO ID RECORC FOR OUTPUTTED LARSYS DATA. 

CONTINUE 

OUTUNT UNIT NUMBER FOR OUTFIT TAPE. 

REC NINBER Of RECORDS IN ONE INPLT FILE. 

RCLL ROLL PARAMETER FCR CNE LARSYS LINE. 

TAPMOC TAPE POCE OF INPUT TA^-E. 

TIPE TINE CF JOB IN SECONDS. 

TIMEl START TIPE OF JOB IN PI LLESECCNO S. 

T1ME2 STOP TIME OF JOB IN PILLESECCNOS. 

TRK7 CONSTANT EQUAL TO 'TTRK*. 

ZERO CONSTANT BYTE EQUAL TO CQ HENIOECIPAL. 


CONTINUE 
>ER CHANNE 


OUTCAT 

OUTIO 


REC 

RCLL 

TAPMOC 

TIPE 

TIMEl 

T1ME2 

TRK7 

ZERO 


C4t****4***<ic«*«««***«*«******«*****«****4«*«4t»)|i*«4*«*««4**«*«««4*4«***4i** 
«««««*«•«•««*«««««*«*««««**•«•*«««««** 44 *** 

E SET .INITIAL CONCITICNS 


TIMEl » TIPERfXl 
COUNT * MAXLC 
CALL MOVBYTfZERO. 


LAROATt 


COUMI 


TRANSFER DATA FILES TC DISK 


C** 4* 4 4 ***************************************** 4********** 4 «* 4* 4 *4«***« 

I MOUNT INPUT TAPE ANO TRANSFER DATA FILES TC DISK 

^****44*********4**4«*********4*4i********»*4****4***4*4**44***4*444«***« 

NRITELUt SC50 

9050 FORMAT?* 1505? ENTER INPLT TAPE AND FILE NUMBERS (I5.IAI*) 
REACI1S»90I0) INTAPt INFIL 
9060 FORMATIISt lA) 

WRITEI 16t9C70 I INTAP* INFIL 

9070 FORMATI* IS07C INPUT TAPE IS**I5.* AND INPUT FILE IS'.U) 


FILE NUMBERS (I5.IAI*) 


>073 ENTER NUMBER CF INPUT FILES 
fAPE ( I3*1X.AA)**/* 

TAPE MOOESC 800 MEANS 9TPACK I 


MODE CF*,/, 


1600 BPI*./ 
MODE 3 800 


Spi*i 


vucu runriAi 1 17, i<t I 

WRITEI 16*9C70 I INTAP, INFIL 

9070 FORMATI* IS07C INPUT TAPE IS*,I5,* AND INPUT FILE IS*,U) 

WRITEI 16,9C73I 

9073 fORMATI* H073 ENTER NUMBER CF INPUT FILES ANO MODE CF*,/, 

1 * TAPE II3,1X.AA)*,/, 

2 * FOR TAPE MOOESC 800 MEANS 9TRACK 690 BPI**/, 

3 * 1600 MEANS 9TRACK 1600 BPI*,/, 

A * TTRK MEANS 7TRACK MODE 3 800 BPI* I 

PEAL'1 15,90161 NCHAN, TAPNCD 

9076 FORMATI I3« IX, A6I 

WRITEI 16, 9C75 : NCHAN, TAPPCO 

9075 FORMATI* IS07* THERE ARE*, 13,* INPUT FILES AND INPUT TAPE*,/, 

1 * MODE IS *,A6,* ISFCSCM*) 

IF (NCHAN .LE. MAXCHN) GC TC 79 
WRITEI6,9C17) NCHAN, MAXCHN 
WRITEI 16, 9C77 INCHAN, MAXCHN 

9077 FORMAT?* E9C77 *,I6,* CHANNELS REQUESTED FCR PROCESSING BUT*,/, 

1 • SPCSCN IS SET UP TO HANCLE CNLV*,I6,* CHANNELS',/, 

2 • PROCESSING ABCRTEO (SPCSCN)* I 

STOP, 

79 CONTINUE 


INPUT TAPE* ,/, 


9077 


SET HEXIOECIPAL 


FLAG ACCORDING TC INPUT 


MODE 


HEXFLG » HEXFf 

IF (TAPMOC .EC. TRK7) HEXFLG * HEX3F 
CALL MOUNT! IN TAP. INUNT, *R0*, TAPPCO) 
IF (INFIL .EQ. 1) GO TO 9C 
CO 80 J = 2, INFIL 
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_ CALL TOfFPlINUNTI 
80 CONItNUE 
90 CONTINUE 


SET UP VARIAeLEE TO TRANSFER FIRST SPECSCAN FILE TO 
DISK FILE 'SPEC FILEl* 

CFILE ■ FILNT 
IDFLG > .FALSE. 

CO 150 J « 1. NCHAN 
REWIND CFIIE 
REC • C 

100 CONTINUE 

READ IN ONE PECCRO FROH TAPE 

COUNT « INBCKT 
REC * REC « 1 

CALL TOFRDIINUNT. COUNT, ERR, INBUFl 

CHECK FOR ENC OF TAPE FILE 

IF I ERR .EO. 1) GO TO 120 
IF (REC .6T. 1) GO TO 105 

SET FLAG IF WE HAVE AN ID RECORD 

CFKIt • C . , 

CALL HOVBYTIINBUF, 0. 1, CHKIO, 3, 1, 11 
IF (CHKID.EO. HEXFLG) IDFLG « .TRUE. 

105 CONTINUE 

IF (ERR .EC. 0) GO TO 110 
WRITEI6, 91001 ERR, REC, J 


91C0 


WRITE! U,910C1ERR, REC, J 

FCRNATI* E91C0 TOPRO ERR«*,I5,* CK RECORD*, 15,* OF FILE*, 


CONTINUE 


(SPCSCM'I 


WRITE RECORD TO DISK 


WRITEI0FILE.911CI INBUF 
9110 FORMAT! )0(1C0A1I,8A11 

GO TC 1(0 
120 CONTINUE 


TRANSFER COMPLETE 


9120 


WRITE! C,9120jJ 

F0RMAT(**I<120 FILE*,I5,» TRANSFERRED TC DISK (SPCSCNI* 1 


: ENC OF TAPE FILE NAS READ SET UP VARIABLES TO 

! TRANSFER NEXT SPECSCAN FILE TC DISK 

CFILE » OFILE 1 
ISO CONTINUE 

CETACH INPUT TAPE DRIVE SINCE ALL INPUT NOW IS CN DISK 

CALL CPFUNC!7,*0ET 182*.CERR1 
IF (CERR .EC. 01 CO TC 17C 
MRITE(6,91(01 CERR 

9160 F0RMAtM*E<16^ IIrCR *,13,* REIURNEO FROM CPFUKC (SPCSCNI*! 
170 CONTINUE 


SET UP OUTPUT LARSVS 10 RECORD AND FOUNT OUTPUT TAPE 


INTRAN IV 
: SPCSCk 
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t081 


t082 


1083 

$384 

1085 

1086 


1087 

1088 

1089 

1090 


! 81 i 


1093 

1096 


1095 

1096 

1097 

1098 


1099 


!ill 

I1C5 

)I06 


1107 


)108 

1109 


1110 

till 


200 CONTINUE 

I INITIALIZE LARSVS IC RECORC 

[ CALL NOVOVIIZEROt Ct C» OLTIOt 0* It 8001 

\ GET OUTPUT TAPE ANC FILE NLRBERS 
HRiTEf Ut9200l 

9200 FORHATI* I <200 ENTER CLTPUT TAPE ANC FILE 
1 *IF YOU ENTER t C lORITE VILL ~ 

REAOtlS.92101 OUTIO(11»OUTIO(2I 
9210 FORFATI I5t 16) 


__ NUMBERS 
SELECT 


l!5jI6)*,< 


FOB YOU*) 


GET RUN NUMBER 


hRITEiUt9230) 

FORMAT (• I <230 ENTER RUN NUMBER TIE)') 
REA0tlSt92<0) 0UTI0f3) 

FORMAT! 18) 


9230 
9260 

STCRE CONTINUATION CODE ANC NUMBER CF CHANNELS 


OUT 1016) 
CUT 101 5) 


0 

NCHAN 


C 

C 

C 


CALCULATE NUMBER OF SAMPLES PER CHARNEL PER LINE 

CALL SPCSAMNSAMP* FILNT* lOFLG) 

0UTIC161 * NSAMP 

GET FLIGHT LINE 


HRITEl IE. 9250) 

5250 FORMAT!* I <25C ENTER FLIGHT LIRE (6A6)*) 

REA0115t92EC) CUTI017)* OUTIOIE)* CLTID19)* OUTIO(IO) 
9280 FORMAT! ^6) 


GLT CATE OAT/! utS TAKEN 


9270 

9275 


WKITEI Ut927C) 

FORHATI* I<270 ENTER MONTH 


DAY, YEAR DATA HAS TAKEN 1316)*) 


ruKnnii' citicK runin* uai* tcak 

REAC115f9275)0UTI0111)f 0UTlC112)t CUTIC113) 
FORHATI 316 I 

SENSOR PLATFORM* 


; GET TINE DATA MAS TAKENt ALTITUDE OF 
\ ANC GROUND HEADING 

NRITE(16f9277) 

9277 FORMAT?* IS27T 


ENTER TIME TAKEN, ALT. OF PLATFCRM* ANC GPOUNO*t/« 
1 • HEALING lA6tIlCfI5)M 

READ115»92T8) OUTIOllA), OUT10115)t 0UTID116) 

9278 FORHATI A6t no tI5) 


GET DATE DATA NAS REFORMATTED 

CALL GTCATIfOUTIDI17) ) 

CALCULATE NUMBER OF LINES CF DATA THAT 
LARSYS RUN —FIRST DETERMINE NUMBER CF 

0UTI012C) 


klLL BE PRODUCED IN THE 
INPUT DATA RECCRCS 


C 

C 

C 


C 

C 


IF llOFLG) Oullo72i) » 0UTI0I2C) - 1 

THEN DETERMINE HOW MANY INPUT RECCRCS MAKE UF CNE OUTPUT LINE 

ONEOUT « INSAMP ♦ IMAXIN -1)) / MARIN 
0UTI012C) « OUTI0120) / ONEOUT 

GET FRAME CENTER LATITUDE AND LONGITUDE 
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9270 FoUatI* ENTER FRAPE CENTER LAT C LCAG I2FT.2I*I 

KEA0I1S.92C1I 0UTI0I21), 0UTI0f22l 
.9281 FORPATI2F7.2I 

\ PRINT OUT AS5ENELE0 ID VALLES FCR RCVIEL 

NRITEI ICtOieOI IQtTI0IJJI.JJ«lt22) 

9280 FORMATI//.* 19280 OUTPLT 10 IS SET UP ASO . ISPCSCM*t/f 

1 • TAFE NC...*iI5»* FILE NC. • I4f/. 

2 • RUN NO....*tI^t* CCNT. CCOE. • •I2#/. 

3 • NO. OF CHANNELS.. .Slit* AC. CF SAPPLES.. . ' t IlC*/* 

A • FLIGHT LINE...*tAAA|/t 

5 • NOATH. DAY. TEAR. ..SI3tl3tI5t/t 

E • TIPE DATA TAKEN. ..* .AA «/t 

7 • ALIITLDE AND GR. HEAOI AG. . . * « I 1 Ct 15*/ t 

E • OATF REF0RPATTE0...S3AAt/. . 

9 • NUPBER OF CATA LINES.. .'tl St/t 

A • LAT £ LONG... *t2F7.2l 


me 


NRITE(Ut929C) 


9290 

FORMAT! • 19290 IS 



READ! 15.92951 ANS 

)ill 

9295 

FORMAT! AA) 

)122 


IF (ANSW .60. NO) GO 


C MCLNT OUTPUT TAPE AND 

H23 


CALL 10RITE(OLTIDt 0 


OUTPLT ID OKAY 


lYES OR AC)*) 


PROCESS CATA —LINE BY LINE 


j CFECK IF OUTPUT NEEDS TC BE FLIPPED LEFT TO BIGHT FRCP INPUT DATA 

9296 FCRPATC** IS29< SHOULD THE OUTPUT CATA BE FLIPPED RIGHT TC LEFT • ,/ 
1 .* ANSWER YES OR NO (SPCSCN)*) 

REACI15t92SS) ANSN 
FLIP « .TRUE. 

IF lANSU .EC. NO) FLIP .FALSE. 

IF tFLIP) VRITE16. 92S7) 

IF IFLIP) VRITEIiet9297) 

.9297 FORMATI//t* 1*297 DATA WILL BE FLIPPED LEFT TC RIGHT (SPCSCN)* ) 

; PREPARE FOR LINE PROCESSING 

CFILE =■ FILNT 
CO 296 J » It NCFAN 
REWIND CFILE 
CFILE B CFILE ♦ 1 
296 CONTINUE 

LCNT - A « (NSAMP * NCHAN) 

NCAT > (NSAMP - 6) * 2 
NCAL « 12 

CHECK FOR OUTPUT GREATER THAN SIZE CF UlTPUT BUFFER 

IF (LCNT .CT. MAXLC) GD TO A20 
NLINE B 0UTI0(2?) 

\ SKIP SPECSPAN It RECORDS IF THEY ARE PRESENT 


IF (.NOT. ICFLG) GO TCi 3CC 
CFILE « FILNT 
CO 298 J B If NCHAN 

REAO(CF lLE.9300fENCB40C) INCUF 
OFILE « OFlLE ♦ I 
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298 CONTINUE 
300 CONTINUE 

PKCCESS DATA LIAE BY LIKE 
CC 390 CURL IN ■ ItNLINE 
READ IN ALL THE CHANNELS FGK ONE LIKE 


CISP • C 
CFIlE « FUNT 
CO 320 CURCHN - 1, 


NCI- AN 


OeiAIN ALL TFE INPUT RECORDS THAT MAKE LP CK£ OUTPUT CHAKKFL FOR 
ONE LINE 


93u0 


DO 31C J « If OKEOLT 

REAO(QFlLE,93COtENO«90C) INBLF 
FCRHATt ICItOCAllfSAll 


TRANSFER THE DATA VALUES TL INPUT DATA LINE EUFFER 

SKIP FIRST 8 BYTES SINCE TI-EY ARE X AND Y INTERVAL 
OFFSETS AND NOT DATA. 

MCOUKT « HAXIN * 2 

If IJ ,EQ. OKEOLT) MCOtNl * KOAT - IMJ-1) * MAXIM ♦ 2) 
CALL HOVBYTI INBLF* €. 1* INCAT. CISP. 1* MCCUNT) 

OISP « OISP « HCCUNT 
10 CONTINUE 

SET calibration BYTES FOR THIS CHANNEL 

CALL PO\BYT(Z£RCf C* C* INDAT. CISP, 1* NCAL ) 

DISP » CISP ♦ NCAL 

CC TO NEXT CHANNEL 


DFILE 

CONTINUE 


DFILE ♦ 1 


9350 


REPACK THE DATA FROM A S-BIT FIELD IN 2 BYTES IMU A 
8-BIT FIELD IN 1 BYTE 

CALL SPCOATTNSAMP, NCHAN) 

SET UP LINE AND ROLL PARAMETER 

LARCAMl) > CURLIN 
LARCATI2) « ROLL 

TRANSFER DATA TC OUTPUT ARRAv, 0ELE1ING NCN-CATA BYTES 
CALL HOVBYTfOUTDAT, 1, 2. LAROATt A, 1. LCNT-A) 

OUTPUT ONE LINE 

COUNT • LCKT 

CALL TOPNRtOUTUNT* COLNT, ERR* LAPCAT) 

IF (ERR .EC. 0) CO TO 360 

WRITE! (*9350 ERR* CURLIN 
WRITH 1(*S35C) ERR* CURLIN 

50 FORMAT! » E9350 TOPHR ERR»'#I5** V.RIT1NC CUTPUT LINE* *16* 

1 • (SPtSCN)*) 

50 CONTINUE 

IF !MOC(CURLIN*100) .EC. C) hR ItE ! 6 .9360) CURLIN 


9360 


390 CONTINUE 
CO TO <90 


IF !MOO!CURLIN.100) .EC. C) WR HE ! 16.9360) CURLIN 
FORMAT!* IS360 ♦,I7,*LINES PROCESSEC (SPCSCN)*) 
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OlfpUl'’YI LEPCTH so LINE COUNT IS NROAC IN 


CONTINUE 

MITff EtV<COINUNEt CURLIN 
NHITEf lE.9<0CINLINEi CURLIN 

FOKPATf* E<AOr. THERE »>AS A PREPATURE ENC-CF-FILE GN CNE CF THE* 

1 ,* INPUT FILE*«/f 

2 ^ SC THE LINE COUNT IN THE 10 0 F*,I 5 »* IS PROeABlV** 

I • hRCNC*./. 

<1 • PROBABLE LINE COUNT IS l*tIS»« • 1) BUT CHECK *t 

i *CUTFUT FILE ANC EDIT THE OUTFUT *•/, 

i • ID IF NECESSARY (SFCSCM't 

€0 TO 4S0 


OUTPUT DATA LINES TOO BIG FOR CUTPUl BUFFER 


A20 CONTINUE 

kRITEIA. 9A20I LCNTt 
l«RiTEI16»9A|0) LCNT. 


I«RITEI16»9A20) LCNT. PAXLC 

9A20 FORMATI/.* E9A20 INPUT DATA NEE0S*tI7»* BYTES CF BUFFER SP/CE BUT* 

1 1* OUTPUT BUFFER IS SET UP*,/. 

2 ' TO HANDLE ONLY’,17,* BYTES PROCESSING**/. 


iw rni.i/i.r uni.T'.ii.' oiic 
7EKPIN/TE0 ABNORMALLY ISPCSCN 


ENC THE JOB 


A50 CONTINUE 


ENC-OF-FILE 

END-OF-TAPE 


HARK AT THE END OF IHE CUTPUT RUN ANC 
RECORD AFTER THAT 


CALL T0PEFI0U1UNTI 
CALL EOTIOLTIC* GUI 


CALL EOTIOLTIC* CUTUNT. ERRI 
PRINT OUT TIPINC INFCRMATICN 
TIME2 » TIPERfXl 

TIME - FL0ATITIME2 -TIMEIl / lOCO.O 
NRItEf 6.9AECI TIME 
WRITE! 16.9A60 I TIME 

9A60 FORMATI/.* I9>60 CPU TIME FOR THl ! JOB WAS**F10.3t* 
1 *ISPCSCN1*I 
STOP 
ENC 


CMftAft IV 6 LEVEL 
E SPCSCM 
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61 

SPCSCN 


PURDUE / 


DATE • 60121 
IAR$ 3031 


10/53/1 


:S1P 


COPEON BLOCK /SPCCCP / PAP SIZE RC41 
LOCATICK SYMBOL LOCATION SYMBOL LOCATICN 

0 OUTOAT AE20 FLIP 9CA0 


SYPeCL 


ir;c< 


KBCL 

MEB 

PAC 

COAT 


LOCATICN 


{l! 


SYMBOL 

MOVBYT 

CPFUNC 

TOPWR 


SUBFROCRAHS CALLED 

LOCATION SYMBOL 

I7B IBCCP« 

18C SPCSAP 

lAO TOPEF 


LOCATICN 

17C 

190 

lAA 


SVPBCL 

PCUNT 

GTOATE 

ECT 


LOC( 


MBCL 

ANK 

XFF 

XLC 

PPCC 


1?“’ 


SCALAR PAP 

LOCATICN SYMBOL LOCATION SYMBOL LOCATICN SVPBCL LOC. 


E9 

FlUNT 

2E6 

F2UNT 

2 EC 

F3UNT 

!§ 

INBCNT 

2FC 

INUNT 

3C0 

PAXO-N 


NO 

310 

CUTUM 

314 

TRK7 


COUNT 

324 

INTAP 

328 

INFIL 

19 

HEXFLG 

338 

J 

33C 

Cf ILE 

'£§ 

ERR 

34C 

CHKlD 

350 

CERR 


JJ 

360 

ANSh 

364 

ERROR 

NC a 

374 

MINE 

376 

CURLIN 

BA 

96 

MCOUN r 

380 

TIME2 

3 6C 

TIME 


MBCL 

TIC 


TSil! 

9075 

916' 

925i 

92 7 

9295 

9360 


•CPTICI 

•CPTIC 

• ST AT 

• STAT 





^RR^Y PAP 





LOCATION 


SYMBOL 

LOCATION 

SYPEOl 

LCICATICA 

SYPBFL 

LOC. 

398 


LAROAT 

6B8 

INBLF 

2CC8 






fCK^AT STATEMEM 

MAP 




LOCATION 


SYMBOL 

LOCATION 

SYMRCt 

LPCATICK 

SY^eCL 

LOC. 

31C7 


9060 

31FB 

9C7C 

32CI 

9073 

3. 

3313 


9077 

3166 

910C 

31E3 

9110 

3' 

346E 


9200 

3AA3 

921C 

35Ct 

92 30 

3* 

3536 


S26C 

3558 

92?C 

355L 

9275 

3' 

35EF 


9279 

35F7 

928 1 

3628 

9280 


3766 


9296 

37HA 

9297 

3819 

93" C 

31 

369A 


940C 

3EC3 

9A2C 

39B6 

9A60 

3. 

NS IN EFFi 

= CT* 

ID tEBCCIC . SOURCE. NOL 1ST 

, DECK, NCLOAC. MAP 



NS IN EFFI 

:CT* 


SPCSCN , LINECN 

t = 

75 






SOURCE STATEMENTS 
NO DIAGNOSTICS GENERATED 


199,PH0gHAP :.ize = 


18339 




ORTRAN IV G LEVEL 21 
E SRCOAT 

C SPCCAT 


62 

SPCORT 


80121 


PUROUE / LARS 3031 


^ SPCCAT TRAKSLITE 9-BIT FIELDS TO C-Bl T FIELDS 

C URITTEA 02/15/79 BY CATHERINE KOZLCVSRI FCR SRCT FV79 CONTRACT 

C REXISEC OA/03/79 BY CATHERINE KQZLCWSKI 

^**««« 4 «*«******«>»« ********»**•««*«***«««*« «**** 4 *«« «*•«««««« 4 4 * 4 * 4 ****« 

C 

SUBRCUTINE SPCOAT (NSAPP. NCh^. 1 
IMPLICIT INTECER 4 A lA-ZI 

C SLBROOTINE TC TPANLATE 9-BtT FIELDS INTC B-61T FIELDS. 

C THE INPUT IS AN ARRAY OF DATA VALUES IN SPECSCAN TAPE 

C FORMAT. 

C EACH DATA VALUE IS ASSUPEO TO BE CNE 3-BIT FIELD IN 

C AN 8-BIT BYTE AND ONE 6-BIT FIELD IN AN 8-BIT BYTE* 

C THESE TWO FIELDS REPRESENTING CNE DATA VALUE RANGING 

C FRDK 0 TO 511. 

C THE OUTPUT IS AN ARRAY CF DATA VALUES WITH EACH 2 BYTES 

C REPRESENTING ONE B-BIT DATA VALUE ITHE FIRST BYTE IS SET TO ZERO 

C ANC THE SECOND BYTE CONTAINS THE DATA I. CUTPL'T VALUES RANGE 

C BETWEEN 0 ANC 255. 

C 

C THIS ROUTINE CALLS NO SUBROUTINES. 

^•«44>4*«*»«««4««44*4*4**«4>«4****4*4>**444*»*44444 4444444*44«4444«444«**44 

C 

c 

COPMON /SPCCON/ 

C INTECER « 2 VARIABLES 

1 INOAT, CUTOAT* 

C LOGICAL * 1 VARIABLES 

2 FLIP 


INTEGER * 2 INOAT I lOOCO) , OUTOATf 1 ?nQC I 
LOGICAL P 1 FLIP 

C444P444444444444444444444444444444444444444444444444444444444444444444, 

€ LOCAL VARIABLES 

C 

INTEGER * A NCAL/6/ 

^4 444 44*444444*44*4 44«4*44i444i44*****44«444i4 444‘44 4444444 44 444 444* 44«**««4 

C LOCAL VARIABLE DESCRIPTIONS 


CCHAN 

CSAMP 

CURIN 

CUROUT 

OISP 


HEIT 

LBIT 

INCAT 

NCAL 

NCHAN 

NCATA 

NSAHP 


CHANNEL CURRENTLY IN PROCESSING. 

SAMPLE CURRENTLY IN PROCESSING. 

DISPLACEPENT INTO INPUT BUFFER FCR CURRENT SAPPLE. 
OISPLACEPENT INTO OUTPUT BUFFER FCR CURRENT SAPPLE. 
DISPLACEPENT INTC LINE BUFFER. FCINTS TC BYTE JUST 
BEFORE FIRST BYTE OF A CHANNEL. 

FLAG INDICATING WHETHER INPUT DATA NEEDS TC BE FLIPPED 
LEFT TO RIGHT ON OUTPUT. IF FLIP IS SET, THE DATA 
SHOULD BE FLIPPED. 

CCNTINUE 

WORK SPACE FOR 3 MOST SIGNIFICANT BITS CF INPUT DATA. 

WORK SPACE FOR 6 LEAST SIGNIFICANT BITS CF INPUT CATA VALUES- 
LINE BUFFER WITH DATA IN INPUT FORMAT. 

NUMBER OE CALIBRATION VAUES PER CHANNL. 

NUMBER OF CHANNELS CF DATA PER LINE. 

NUMBER OF DATA SAMPLES PER CHANNEL. 

NUMBER 01 SAPPLES PER CHANNEL. 


IRTftAM IV 6 LEVEL 
: SPCOAT 


21 


63 

SPCDAT 


DATE > 80121 
PURDUE / LARS 3031 


13/A9/0 


QUIOAT LINE CUFEER LlTH OTA IN OUTPLT FCRMAT. 

TEMP WORK SPACE USED IN CALCULATICN OF MBIT AND LB IT. 


1008 

1009 

1010 
1011 


1012 

1013 

lOlA 


1015 

1016 
•017 


tOlB 

>019 


020 

021 

022 

023 

024 




C 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


s 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


SET NUMBER OF NON-CALIBRATION OATA VALUES 
NOATA « NSAMP - NCAL 

PROCESS the input ARRAY CHANNEL BY CHANNEL 
DO 400 CCHAN « It NCHAN 

CALCULATE POINTER TO BEGINNING OF CURRENT CHANNEL 
OISP « (CCFAN - 11 « NSAMP 
PROCESS THE CAT I CHANNEL SAMPLE BY SAMPLE 
CO 300 CSANP s It NOATA 

CALCULATE PONTERS TO CURRENT INPUT CATA VALUE ANC 
ITS CORRiSPONOING OUTPUT DATA VALUE 

NOTE THAT OUTPUT DATA MILL BE REVERSED LEFT TC 
RIGHT FROM INPUT DATA IF FLIP FLAG IS SET 

CURIN « OISP * CSAMP 
CUROLT « CURIN 

IF IFLIFI CUROUT » OISP ♦ CNOATA - CSAMP ♦ 


II 


PUT THE 2-BIT INPUT FIELD INTO THE 
TEMPORARY OUTPUT 9-BIT FIELC 

TEMP « INDATfCURINI 
HBIT > TEMP / 256 
♦BIT » HBIT * 64 

PUT THE 6-BIT INPUT FIELD INTO THE 
TEMPORARY OUTPUT 9-BlT FIELC 

LBIT * NOOfTENP. 64) 

TEMP = HBIT ♦ LBIT 


3 MOST SIGNIFICANT BITS OF 


6 LEAST SIGNIFICANT BITS CF THE 


IN THE CUTPUT ARRAY. 

INTO AN 8-BIT FIELD BY ORCPPING 


STORE THE TEMPORARY FIELD 
TRANSLATE THE 9-BIT FIELD 
THE RIGHT-MOST EIT. 

SINCE THE INPUT IS A NEGATIVE IMAGE. THE CUTPUT IMAGE 
IS INVERTED TO MAKE A POSITIVE IMAGE. ITHAT IS. A TEMPCRARY 
FIELC OF 0 IS OUTPUTTED AS 255 AND A TEMPORARY FIELD 
OF 255 IS OUTPUTTED AS C. I 


OUTCATLCUROUTI 
300 CONTINUE 
4C0 CONTINUE 
RETURN 
END 


» 255 - fTEMF / 2) 


OMTRAN IV 6 LEVEL 
E SPCSAf* 


2t 


64 

SPCSAP 


3001 

3002 


3003 

3004 
300$ 
3006 
3037 


3008 


PtROUE / 


DATE * 80121 
LARS 3031 


13/49/0 


SPCSAH 


SPCSAM CLACULATE NUMBER CF SAMPLES IN LARSVS OUTPUT 
MRtTTEN 04/03/7S BY CATHERINE K0ZL06SKI FOR SRST FV79 CCKTRACT 
«M*«4***4««««»««««**«*9**********»*«*«9««««««*4«*«««««*«««««4«4*4**«*4t 
SU6RCUTINE SPCSANINSAMPt UNIT, lOFLG) 

^ IMPLICIT INTECER * 4 lA • 2) 

C«**P44P**«****P4«*«*****«*«*««*«**«***«*<»*4«**«4*««#«««««#*«««* <«««*«** 

C . SPCSAM DETERMINES THE NUMBER OF SAMPLES FER LINE FER CM^NNEL 
C IN THE LARSVS FCRMATTED OUTPUT RUN* THE RCITINE READ INPUT 
C RECORCS OFF CIS* UNTIL AN INPUT RECCRD ENDS VITH ZERCES. 

C T*E ROUTINE THEN COUNTS THE NUMBER CF NCN>ZERO INPUT VALUES THAT 
C HAVE BEEN READ INO RETURNS THAT COUNT PLUS CALI6RATICN BYTES AS THE 

C NUMBER OF SAMPLES PER CHANNEL. THE INPUT TC THE ROUTINE IS THE UNIT 

C NUMBER OF A CISK FILE CONTAINING INPUT DATA PECCRDS AND AN INPUT 
C IC FLAG. THE OUTPUT FROM THE ROUTINE IS THE NUMEER OF SAMPLES. 

C THIS ROUTINE ASSUMES THE DISK FILE FAS ALREADY BEEN SET UP. 

{ THIS ROUINTE CALLES ONLY SYSTEM SUBROUTINES* 

C444)*44«*4*4i*«4>««*4>«4c4(4i*4t4E**4i*«4i4>«4t«*«*44>4i4i«44i««4««4««**«4««**4« «**«***4i 

c 

INTEGER * 4 BLFCNT /5C4/, MAXOAT /500/t NCNOAT /4/, RLIMIT /$/, 

^ 1 ZERO /O/ 

INTEGER « 2 BUFFER 1 5041 

^ LOGICAL « 4 ICFLG 

^ LOGICAL * 1 INBUFflOOfil 

EQUIVALENCE CEUFFER, IN6LFI 

c 

C LOCAL VARIABLE DESCRIPTIONS 

^ ADJUST temporary VARIABLE USED TO ACJUST SAMPLE COUNT TC BE EVENLY 
C DIVISIBLE BY 4. 

C 6UFCNT NUMBER OF ELEMENTS IN INPUT 8LFFER. 

BUFFER INPUT BUFFER IN INTEGER • 2 FORMAT. 

OISPLACENENT INTO INPUT BUFFER. 

FLAG INOICATIANG WHETHER FIRST INPUT RECORD CF THE FILE IS 
AN ID RECORD. IF lOFLG IS SET. FIRST RECCRD IS AN ID. 

INPUT BUFFER IN LOGICAL * 1 FCRMAT. 

LAST CAT A VALUE IN INPUT BUFFER. 

MAXIMUM rUMBER OF DATA VALUES POSSIBLE IN ONE INPUT RECORD. 

CONTINUE 

NCNCAT NUMBER OF NON-DATA VALUES IN INPUT BUFFER. 

NREAO NUMBER OF CONSECUTIVE RECORDS READ WHEN SEARCH FCR ZERC DATA 
VALUES. 

NSAHP NUMBER OF SAMPLES PER CHANNEL THAT WILL BE IN LARSYS CUTPUT. 
PREVAL PREVIOUS DATA VALUE IN INPUT BUFFER. USED TC SEARCH 
BACKWARDS IN INPUT RECORD. 

RLIMIT UPPER LIMIT ON THE NUMBER CF CONSECUTIVE REACS TC PERFORM 
BEFORE TERMINATING SEARCH FOR ZERC DATA VALUES. 

TOTAL RUNNING TOTAL USED TO CALCULATE NUMBER CF SAMPLES. 

UNIT DISK UNIT FROM WHICH TC READ INPUT RECORDS. 

ZERO THE CCNS1ANT ZERC. 


DISP 
ICFLG 

INEUF 

LSTVAL 

MAXCAT 


CP*«*4 4«*4**4m>«*«** 4*«*44«****<» 4i*4i 4i**«4>4««* 44**4 4<*«*444«44* 4**44 4444>***« 
^4444*4**4«*4***444«*4444«44**4**«*«***44444 4*4 44 *4*444444*4 4 4 4* 444**4 4* 


d 


IRTIIAN IV 
: SPCSAF 


)0C9 

tolo 


>911 


)012 

>013 

>014 

>015 

>316 

>017 

>018 


1014 

>020 


>021 

>022 

>023 

>024 

>025 


>026 

>327 

>028 

>029 

>030 

>031 

:sii 

>034 

>035 

>036 


>037 

>038 

>039 

>040 

>041 

>042 


6 lEVEL 21 


63 

SPCSflH 


PURCUE / 


DATE « 80121 
LARS 3031 


13/49/0 


INITIALIZE VARIABLES 

: 8 


NREAI 

^«**4«******4***««««***4>*********i|i***«***«*«««**«*«««««**4*«4«4* <««*«**« 

READ CATA FRCN FIRST DATA FILE UNTIL DATA VALUES EQUAL ZERO 

;«««*«#««9«*«««*«**«*4i4>**«*«*«*«««*«*«**«i|i*4«**«4*4*«*»**«««*«4* <**«*«** 

REWIND UNIT 

SKIP ID RECORD IF THERE IS ONE 


91 


tl8 


IF CIOFLGI RE/OIUNIT.SIOC INBUF 
FORKATI lOnOOAl 1 1 BAl I 
CONT INUE 

REACIUNITtSlOOl INBUF 
LSTVAL a BIFFERIRUFCNTI 
TOTAL « TOTAL * MAXOAT 
NREAO « NREAD ♦ 1 


C 

c 

c 

c 

c 

c 

c 

c 


STOP READING RECORDS IF ZERO VALUES FOUND OR RLIPIT RECCROS HAVE 
BEEN READ IN A FOW 

IF CILSTVAL *NE. ZEROI .AND. INREAO *LT. RLIPITII GC TC I2C 


IF ILSTVAL .EC. ZEROI 60 TO 200 

THERE WERE NC ZERO CATA VALUES IN RLII4IT COXtSECUTIVE CATA RECORDS 
ASSUME 1 SPfCSCAN RECORD EQUALS I LARSVS CHANNEL PER LIKE 

TOTAL a MA>DAT 
WRITEI6. 9150) 

WRITE! lEt91S0l 

9150 FORMAT!/.* IT IS ASSURED THAT 1 SPECSCAN RECCRC EQUALS* » 

1 • 1 CHANNEL PER LINE OF LARSTS DATA fSPCSAI'I'l 

GO TO 290 

t COUNT NUMBER OF ZERO DATA VALUES ANC DELETE FRCP TOTAL CCUN'T 

2C0 CONTINUE 

CISP a eUFCNT 
PREVAL a ZERO 

220 CONTINUE 

IF flOISP .LE. NONCAT) .OR. IPREVAL .KE. ZERO) CC TC 240 
PREVAL a BUFFER(OISP) 

IF IPREVAL .EQ. ZERO TOTAL » TCTAL - 1 
OISP a CISP - 1 
60 TC 22C 

240 CONTINUE 

IF lOISP .CT. NONDAT) GO TO 27C 


C 

C 

C 


THERE ARE INPUT RECORDS WITH ONLY ZEROS IN THER 

INCLUDE SOME OF THESE ZERO DATA VALLES IN THE OUTPUT SC 
THERE WILL BE PROPER CORRESPONDENCE BETHEEN INPUT 
RECORDS AND LINES OF OUTPUT. 


e-240 


29 


WRITEI6* 9240) 

WRITE! K. 9240) 

FORMAT! INPUT INCLUDES DATA RECORCS THAT CCKTAIN* , 
I * ONLY ZEROS ISPCSAR)*) 

TOTAL a TOTAL ♦ 10 
CONTINUE 
CONTINIE 


C SET UP OUTPUT SAI^PLE COUNT 



i 

i 





1 
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68 

X1CCR4 


DATE « 8C121 
PURDUE / LARS 3031 


10/53/1 


X10CM4 CONTAIAS ALL VARIABLES NEEDED 1C REFORPAT EXCTECH HCOEL 

IOC DATA. A SEPARATE COPKCN RLOCK LAS SET UP TC KEEP THESE 
VARIABLES SEPARATE FRCP 1HF VARIABLES USED FCR REFORPATTING 
OF EX01ECH PODEL 2CC DATA. 


WRITTEN 03/15/79 RV CATHERINE K02L0VSKI 
REVISED 07/19/79 BY CATHERINE K02L0VSKI 


BLOCK DATA 

IMPLICIT INTEGER • A IA-21 


COPPON /XICCP^/ 
REAL * * VARIAELES 


1 CLR 

2 OKA 

3 6AI 
A CSX 

Integer 


CLRAFTIlElf 
OKAFTf 121* 
6AlN(9lf . 
GSXlOf 12lt 


CLRCALIAtZli 

EXPARBf Alt 

iNSUNTt 

MAXeNDt 

NUPBERf 

PTRRt 


variables 


TIMAFT* 

DLRBEFIl2)t 
DRBEFf 12lt 
GRAFTI 12 It 
RRXIOI 121t 
THRXlOt 


0RCALI.At2)t 

EXP ARC! Alt 

LNUPBt 

MAXGANt 

NUMINSt 

PSCONVt 


c rini%t r t 

LCGICAL * 1 VARIABLES 
C BNDSKPI121t OBGCINt 

C CBGSFWt NEWCINt 


IlMBEFt 

CLRXI0(121t 

CRXlCfl2lt 

GRBEFI121t 

IHSXlOt 

CTYPEt 

FNUMBt 

NDANCt 

FRESEKf 

RESSHEt 

CBGSBK, 

SKPOBS 


0LSX10I121 t 
OSXlOf 121t 
6RX10tl21t 

X10CNC(12t21t 

EXPARACAlt 

INSTNAIAlt 

NGAIN. 

PRT5I21t 


OBGSCHi 


real * 8 

.REAL * A 

1 CRAFTt 
3 GAIN. 

A GSXlOt 
5 RRX lOt 

INTEGER « A 

1 

2 EXPARBt 

3 INSUNTt 
A MAXBNOt 
5 NUMFERt 
t PTRRt 

LOGICAL « I 

1 

2 CBGSFWt 
EQUIVALENCE 


TINAFTt 

OLRAFTt 

DLRBEF t 

ORBEFt 

GRAFT. 

REXPAAIAlt 

THRXlOt 

DLRCALt 

ORCAL. 

EXP ARC f 
LNUPB. 
MAXGAN f 
NUMINSt 
PSCONVt 

BND SKP t 
OBGCINt 
NEWCINt 


TIMBEF 


WillV 

GRBEFt 

REXPABIAlt 

THSXlOt 


CTYPEt 

FNUMBt 

NBANC. 
FRE SERf 
RESSHE 


CBGSBKt 

SKPOBS 


OLSXlOt 

OSXlOt 

GRXlOt 

REXPACIAlt 

XIOBNC 


EXPARAt 

INSTNAt 

NGAINt 

PRT5t 


OBGSCI-t 


(EXPARAIllt REXPAAflllt lEXPARBIll tREXPABCll It 
(EXPARCf lit REXPACtll 1 


X10CM< VARIAELES 


1 


ilRTRAN IV 6 LEVEL 21 
r XlOCf'4 


69 

BLK CATA 


\ DATE « 60121 

PURCUE / LARS 3031 


BNCSKP FLAGS DATA VALUES TC SKIP WHEN PROCESSING AN ^ 
OBSERVATION. IMTIALI2EO HERE TO PEAN NC INFCR- 
HATICN KNOWN VET. RESET IN SUBRCUTINE KIORCL PROP 
CALIERAltON SHEET. 

CATA BNOSKP /12 « F/ 


OEGCIN DEBUGGING FLAG. IF SET ON* DEBtG INFCRPATICN IS PRINTED 
FOR SUB POUT INE KIOCIN. 

OEGSBK DEBUGGING FLAG. IF SET ON* OECL6 INFCRPATICN IS PRINTED 
FOR SUBROUTINE K10S6K. 

CBGSCH DEBUGGING FLAG. IF SET ON* OEBLG INFCRPATICN IS PRINTED 
FOR SUBROUTINE KIOSCH. 

*^BGSFW CEBUCGING FLAG. IF SET ON* OEBLG INFORPATION IS PRINTED 
FOR SUBROUTINE XIOSFW. 

OLRCAL CARK LEVEL CALIBRATION INSTRUCT ICNS. INITIALIZED 

HERE TO MEAN NO INFORMATICN KNONN YET. RESET IN THE 
SUBRCUTINE XIORC ACCORDING 1C THE CAL16RATICN SHEETS. 

CATA OLRCAL /f*0/ 

OLRAFT DARK LEVEL REFERENCE FOR REFERENCE STANCARO. USEC ONLY 
FOR CALIBRATION COCE 23. OBTAINED PROP SUBROUTINE 
XICOREF ACCORDING TO CALIBRATION INSTRUCTIONS. 

OLRBEF DARK LEVEL REFERENCE FOR REFERENCE STANDARD. USEC ONLY FOR 
CALIBRATION CODE 23. OBTAINED PROP SUBRCUTINE XIOREF 
ACCORDING TO CALIBRATION INSTRUCTIONS. 

OLRXIO CARK LEVEL REFERENCE FOR REFERENCE STANDARD. OBTAINED PROP 
SUBRCUTINE XIOREF ACCORDING TO CALIBRATICN INSTRUCTIONS. 


CBGSCH 


*^6GSFW 


OLRCAL 


OLRAFT 


OLRBEF 


OLRXIO 


OLSXIO 


CRCAL 


DARK LEVEL REFERENCE FOR SCENE. 
SUBRCUTINE XIOCLM. 


CBTAINEC FRCP 


REFERENCE STANDARD CALIBRATICN INSTFLCTICNS. INITIALIZED 
HERE TO MEAN NO INFCRPATICN KNOWN YET. RESET IN THE 
SUBRCUTINE XIORCL ACCORDING TC THE CALIBRATICN SHEETS. 

CRCAL /6«C/ 


DRAFT REFERENCE STANDARD CATA RESPCNSES. USEC ONLY FCR 

CALIBRATION CODE 23. OBTAINED PROP XIOREF ACCORDING TO 
CALIBRATION INSTRUCTIONS. 

DRBEF REFERENCE STANDARD DATA RESFCNSES. USEC ONLY FCR 

CALIBRATION CODE 23. OBTAINED FROP XIOREF ACCORDING TO 
CALIBRATION INSTRUCTIONS. 

DRXIO REFERENCE STANDARD DATA RESFCNSES. CBTAINEC FRCP 

SUBROUTINE XIOREF ACCORDING TO CALIBRATICN INSTRUCTICNS. 

CSXIO SCENE DATA RESPCNSES. OBTAINED FROP RESPCNSE 

SHEET. (ALSC CONTAINS PRT-S* EXPARA* EXPARB* AND 
EXPARC RESPONSES.! 

OTYPE DATA TYPE FOR SCENE. INITIALIZED HERE TC PEAN NC 

INFORMATION KNOWN YET. RESET IN THE SUBROUTINE XIORRS 
FROM DATA RESPONSE SHEETS. 


CATA DTVPE /• 


3RTRAK IV G LEVEL 
: X10CR4 


21 


70 

BLK DATA 


>013 


)014 


>015 


>016 


C 

C 

C 


C 

{ 


c 

c 


c 

I 

c 


c 

c 

c 


PtROUE 


DATE • 80121 
/ LARS 3031 


10/S3/1 


EXPARA ARRAV CCNTAINING INFORPATIDA FOR AN EXPERIPENTER PARIMETER. 
EXPARAflt CONTAINS THE NUPBER OF 1HE ID ELEPEKT FCR THIS 
EXPERIMENTER PARAMETER* EXFARAI2I CCNTAIN THE KUPBEF OF THE 


DATA REIPCNSE CHANNEL FRQP AHICH THE VALLE SHOULD COPE. 
|XPARA(3) CONTAINS THE PRT-5 CONVERSION TABLE 
THE REP!nSE. EXPARAIAl STORES THE CALCULATED 


VALUE IN REAL FORMAT. EXPARAllI THRCUGH 
13) CBTMNEC FRCM THE SUBROLTINE XICRRS. 

SET I 

NOTEC _ _ _ 

AGRONOMY SHEETS THE AGRCNCPY SHEET VALUE IS USEC. 


EXPARAIAI 

IN SUBROUTINE XICCLM. 

If THE EXPERIMENTER PARAMETER VALUE IS ALSO IN THE 


EXPARB 

EXPARC 

FNUMB 


GAIN 


SAME AS EXPARA. 
SAME AS EXPARA 


NUMBER CF FIRST OBSERVATION TO VHICH THE CURRENT CALIBRATION 
INSTRUCTIONS APPLY. OBTAINED FROM SUBROUTINE XIORCL 
ACCORDING TC CALIBRATICN SHEETS. 


GAINS TO USE IN REFORMATTING DATA, INDEXED BY GAIN CEDE. 
INITIALIZED. HERE. TC MEAN.KO INFCRMATICN KNOUNYET. 


RESET 


SUBROUTINE XIORSS ACCORDING TC RESPONSE SHEETS. 


CATA GAIN /V«C.O/ 


GRAFT 

GR6EF 

GRX13 

GSXIO 

INSTNA 

INSUNT 

DATA 


GAIN USED FOR REFERENCE STANDARD. _ 

CODE 23. OBTAINED FROM SUBROUTINE XIOREF ACCCRCING TO 
CALIBRATION INSTRUCTIONS. 


USED CNLY FCR CALIBRATICN 


GAIN USED FOR REFERENCE STANDARD. 
CODE 23. OBTAINED FROM SUBROUTINE 
CALIBRATION INSTRUCTIONS. 


LSED CNLY FCR CALIBRATICN 
XIOREF ACCCRCING TO 


GAIN USED FOR REFERENCE STANOARC. CETAIKEO FCR THE 
SUBROUTINE XIOREF ACCORDING TO CALIBRATICN INSTRUCTICNS. 

CAIN USED FOR THE SCENE. OBTAINED FROM DATA RESPCNSE 
SHEET. 

THE INSTRUMENT NAME OBTAINED FRCP THE SUBROUTINE XIORSS. 
UNIT NUMBER FOR INSTRUMENT CODE TABLE. 

INSUNT /17/ 


LNUNB NUMBER CF LAST OBSERVATION TC WHICH THE CURRENT CALI ERAT ION 
INSTRUCTIONS APPLY. OBTAINED FROM XIORCL ACCCRCING TO 
CALIBRATION SHEETS. 

MAXBND MAXIMUM NUMBER OF BANOS THE SYSTEM CAN HANDLE. 

CATA MAX8NC / 12/ 

MAXGAN MAXIMUM NUMBER CF GAIN COOES. 

CATA MAXGAN /«./ 

NBAND NUMBER CF BANOS OF DATA RESPONSES. 

NGAIN NUMBER CF GAINS TO USE FOR THIS SET CF DATA. 

NEWCIN FLAG TO INDICATE A NEk SET CF CALIBRATION 

INSTRUC1I0NS. SET IN THE SLBROLTINE XIORCL AND RESET IN 
THE SUBROUTINE XIOREF. 


IRTDAN IV 6 LEVEL 
: X10CF4 
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BLK DATA 


k DATE « 801?1 

PUROUE / LARS 3031 


NtPBER NUMBER <8 DIGITS) CF THE OBSERVATICN CURRENTLY BEING 
REF0RMA1TEC* OBTAINED FOR IhE SUBRCLTIAE XIORFR. 

NURINS NUMBER OF DARK LEVEL OR REFERENCE STANDARD INSTRUCTICNS. 

CATA NUNINS /4/ 


PRESER CONTAINS THE PREFIX 


OBTAlNEf FROM THE SUBROUTINE XICRSS ACCORDING TC THE CATA 
RESPCNSE SHEETS. 

PRT5I2) ARRAY CCNTAINING INFORMATION FOR PRT5 DATA VALUE. 

PRfTl) CONTAINS THE NUMBER CF DATA RESPONSE CHANNELS. 
PRTI2) WNTAINS THE PRT-S CCNVERStCN TABLE. 

P1RR POINTER INTO RESPONSE SHEET FILE. 

PSCONV UNIT NUNBER FOR PRT5 XlOO DISK FILE. 

DATA PSCONV /24/ 


FOR THE ID SERIAL NUMBER 


RESSHE nSK UNIT NUMBER FCR RESPONSE SHEET FILE. 

CATA RESSHE /-S/ 

RRXIO SPECTRAL BIDIRECTIONAL REFLECTANCE FACTOR OF REFERENCE 
STANCARC. OBTAINEC FROM SUBROUTINE XIOTAB FROM 
DISK FILE RFLSTND XlOO. 

SKPOBS FLAG TO INDICATE OBSERVATION SHOULD NOT BE PRCCESSEC. 
SET IN XIOREF. 

THRXIO SOLAR ZENITH ANGLE AT TIME CF OBSERVATION CF REFERENCE 
STANCARC. OBTAINEC FROM SUBROUTINE XIOREF. 

THSXIC SOLAR ZENITH ANGLE AT TIME CF OBSERVATION OF SCENE. 
OBTAINEC FROM SLBRCUTINE XlCREF. 


TIMAFT 


OBTAINEC FROM SLBRCUTINE XlCREF. 

TIME OF REFERENCE. RESPCNSE AFTER THE SCENE RESPCNSE. 
USED ONLY FOR CALIBRATICN CCOE 23. INITIALIZEC 
HERE TO MEAN NO INFCRMATION KNOWN YET. RESET 
IN THE SUBROUTINE XlCREF. 


CATA TIMAFT /C.O/ 


TIMBEF 


REFERENCE RESPCNSE BEFORE THE SCENE 
RESPCNSi. USED ONLY FCR CAIIBR/TICN CCCE 23. 


HERE TO MEAN NO INFORMATION KNOWN YET. 
SUBROUTINE XIOREF. 


BEFORE THE SCENE 
IIBR/TICN CCCE 23. INITIALIZED 
KNOWN YET. RESET IN THE 


CATA TIMBEF /C99.9/ 


XIOBNO LOWER AND UPPER LIMITS FOR EACH BAND CF EXCTECK 

MODEL ICO DATA. OBTAlNEO IN THE SUBROUTINE XIORSS FROM 
THE INSTRUMENT CODE DISK TABLE. 
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Appendix A-4 


Control Card Reference File 


and 


Program Abstracts 


for 


SMOOTMRESULTS 


and 


CHANGEDETECTION 


Programs 


KtVISEO 01/^4/fO 
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LAkSVS CONTROL CA«<US 
SMOOT M» »:SUl TS 

t KEV CONT><OL 

O MOPO(COLvl) MAHAMETtR f UNCI ION 


♦ •SWOOTHMF.SUITS INONtI 


OtPLftCES Or<UUPS Of POINTS 

<*iT*i Trie Dominant Class. 


♦ inpesults 


TAMt(XXA) 
FUt iff ) 
I?IS- 


LOCaTION of jn*»ut .results 

rnCAfLU Urr TAHt XX A 
MLf. F> 

L'St results HLmCEo on DTSiv 
tor Ct ASSlf rPOlNl> IN CO«r)tNT 

Tc"*minal session. 


CFLLSI/t LL*CC 


UFFInE Tht CELL DIMENSIONS* 

LL Ts The NOMtt*' Of LINES. 

Cc I> 1r»t NOMre**' OF CO| OMNS. 
f^iaxiMO** SWr. IS lu X lo. 


♦ OOTKtSOLTS 


TAPE (XXX) 
FiLtIFE ) 
IimITIalIZF 


niS" 


OFSTInaTIoN or h-ESULTS 
PUT On TAPt XXX. 
file FF. 

iNllIALlxt file one of a 
nEw PESulIS 
rPEOOIPtD *Ht.r4 USlN(» A 
u-Ew fAPt). 

EFSULTs WILL oE STci*^EO ON 
C AWSrS 0 1 bn . 


PklOPITY GI.G^«... 


PulomiY o 

WILL NOT 

The (.ell 


NOUPs Gi.Gf..** 

Ht PEPL^'CHU wr*£N 
IS HuDlFItD. 


GWODP NAPtlGl/Pl*Pa/> GPOOP classification POOLS 

P1*PZ«... FOP CALCtLATl.JG 
COhRElI PtCOGNlTION* 'NAME, » 
IS The owoup name and gi is 
ThE GkOOP NoahEP. 


^EIGHTS 


Wlfw^*... ASSIGN WEIGHTS 10 hoCi.S. 

In T-'Is Cwi'.EP. 


OLOCn 


wi/N (xxxAxxxx ) Data FhO’^ wwh xaxxaxxx is-— 

He. ■ -»tSTfJ. 


tiNtlX.y./) I 'A I A FwOM I iNt A TO r WITH 

jNTr-VAl / 

roi (X.Y.Z) DATA FmOM column X T(J Y * 1 TH 

lUTr-Vttl — — 


default 

(NONE) 

(NONE ) 

(ZX2) 

(NONF ) 

SEE CONTPOL CAPO 

DICTIONaBv 

(NONF) 

NO GPOUPING 

EGUAl. weights used 
ENT I PE AREA 


♦ EnO 


( NONE ) 


E* l* OF fONLlIOt-! 


(NONE ) 


74 LARS Program Abstract 

MODULE IDEWTiriCATION 

Module Mamet smosup Function Name: smootoresults 

Purpose : Suparviaor for SMOOraRESULTS 

System/Language t CMS/PORTiuut 

Author : John Cain Date ; 4/1B/80 

Latest Revisor t Date: 


MODULE ABSTRACT 

SMOSUP supervises the SMOOTORESULTS function by celling two subroutines: 
one to read the control cards, euid one to modify the input results file. 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 
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SMOSUP-2 


1« Moduto u»aq» 

CALL 8M0SUP 

No argmont* ara usad in tha call to SMOSUP. This SMOOTORESULTS supervisor 
routine is called from LARSMN. Upon conplation of tha function, control is 
ra turned to IAR8HN. 

2« Internal Description 

8M06UP contains two cononon blocks • CLOCOM and SMOCOM. This supezvisrry 
routine calls 8M0RDR to read in tha functim control cards. After the 
cards have bean interpreted, SMOSUP calls SMOINT which initiates the 
modification of an area. When SMOINT returns control to SMOSUP, the 
supervisor indicates that the function is completed and returns control 
to tARSMN. Subroutines called by SMOSUP: 

SMORDR 

SMOINT 

3. Input Description 
Not ^^lioable. 

4. Output Description 

Standard supervisor information messages. 

5. SupplMtental Information 
Not sg>plicable. 

6. Flowchart 


Not applicable. 
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MODULE IDENTIFICATION 

Module Name; smooom Function Name: smoothresults 

Purpose : r«»iminn fr>i> r.MnOTTO 1 gSlir.Tg 

Sys tem/Language : cms /fortran 

Author: John Cain Date t 4/18/80 

Latest Revisor: Date: 


MODULE ABSTRACT 


I 

'i 


i 


I 





PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47*#C6 
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LARS Program Abstract 

MODOLE IDENTIFICATION 

Module Name: smoint Function Name t smoothresults 


Purpose: Modify the 

input Classificati<») Results File 

System/Language : 

CMS /FORTRAN 

Author : John Cain 

Date : 4/18/80 

Latest Revisor: 

Date: 


MODULE ABSTRACT 


SMOINT reads the input Classification Results Pile "n" lines at a time 
(where n is the number of lines in a 'cell'), reassigns class numbers 
to the dominant class (depending on user-input parameters) and writes 
out a modified Classification Results File. 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
Wrfce»*- 47906 
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SMOINT-2 


1» Module Usage 
SMOIMT 

CALL SMOINT 

There are no arguments to SMOimj tdiich is called by SMOSUP* and returns 
control to SHOSUP when the function terminates. 

2 . Internal Description 

SHOINT reads from the input classification Results File cind modifies each 
record type in the following way: 

RECORD TyPE 1 - Record ^pe 1 is read emd checked to see whether it 
is valid (i.e., not a startup file). Ihe new tape and file numbers 
are written onto the output Classification Results File. 

RECORD TYPE 2 - The number of classification pools is changed to the 
nunber of group names amd the pool pointer and stack arrays are 
clumged accordingly. If no grouping was done, the pools are then 
considered classes in themselves. If there are any weights, then 
they are checked and normalized before written onto the output results 
file. 

RECORD TfPE 3 - The first card of the statistics deck is marked 
indicating that the deck is invalid due to execution of the 
*SMOOTHRESULTS function. This deck is then copied onto the output 
results file. 

RECORD TfPE 4 - unchanged. 

RECORD TfPE 5 - SMOINT checks the area requested by the user to 
be sure it exists. If only part of the area requested exists, then 
the user is given the option to continue or terminate the function. 

If the user continues, the parameters are changed so that they are 
now valid. These are written to the output results RECORD TYPE 5. 

If no BLOCK card was used this record remains unchcmged. 

RECORD TYPE 6 - This record is read into buffers "n" lines at a time, 
(where "n" is the number of lines in a cell) and shifted until only 
the columns requested are considered. SMOINT then calls SMOOTH to 
modify the buffer data. Upon return, SMOINT writes the modified 
lines to the outqiut Classification Results File. 

RECORD TYPES 7, 8: Taese record types are created according to 

the required Classification Results File format. A final record 
TYPE 1 is written, all tapes are detached and control returns to 
SMOSUP. Subroutines called by SMOINT: 


SMOOTH 
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SMOINT- 


3» Input Description 

8M0XNT reads the first 6 record types off of a lAJtSYS Classification 
Results File. It also checks the weights card and can read in a 
corrected version if required. 

4. Outt>ut Description 

Several information messages may be issued. Obey are as follows: 

1 ^) I*** nils is a restaurt file — Run Class! fypoints first. a. 

b) I*** File length is only one record — Try running Listresults. 
q) X*** Area requested only partially within this classification 
area. Do you wish to continue? 

d) I*** Execution terminated by user. Do not consider partial area. 

All ei^t record types axe written to either tape or disk in the LARSYS 
Classification Results File format. 

5. Supplemental Inforamtion 
Kot ^plicable. 

6. Flowchart 


Not applicable. 




MODOLE IDENTIFICATION 
Module Name : s**®®®*^ 


LARS Program Abstract 


Function Name: shootoresults 


Purpose : Read and Interpret c 

System/Language i cms/fortran 

Author : Jbhn Cain 

Latest Revisor : 


Date: 4/18/80 
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SMORDR-2 


1. Module Usage 
8M0RDR 

CALL SMORDR 


there are no argumente to SMORDR. SMORDR interprets the following 
function control cards and takes the indicated actions: 


CONTROL CARD 
INRESULTS 


CELLSIZE 


OUTRESULTS 


PRIORITT 

GROUP 


ACTION 


INRES is set to the correct DSHN for tape or 
disk. If it is a tape that is requested, 

MMTAPE is called to mount and position it. 

CELROW is set equal to the n\amber of lines 

per cell (length) and CEIXTOL is set to the number 

of columns (width) . 

OUTRES is set to the correct DSRN for tape 
or disk. If it is a tape, MMTAPE is called 
to mount and position it. 

The vector PCLASS is filled with user 
defined priority classes. 

GRPNAM and GRPSTK aro conputed by a call to 
GRPSCN. 


WEIGHTS 


These weights are read into PROB^ a REAL*4 
vector of dimension GO. 


BLOCK 


MIXCLASS 


The first 6 words of array BLOCK are used to 
contain this field boundary definition and 
the run number is put in RUNNUM. 

The array MIXDAT is filled with the low and 
high range percentage values to be later 
conpared to each group in each cell. The 
array is filled by a call to READMX. 


2. Internal Description 

SMORDR initializes necessary variables and then uses CTLWRD to interpret 
the keyword on each card. An unexpected end of file for the control card 
input causes a call to ERPRNT and termination of the function. Once tliO 
keyword has been interpreted, SMORDR uses CTLPRM, IVAL, and FVAL to 
'urther interpret the card. If a disk is requested the function TSPACE 
i.s used to calculate the amount of available space on the disk. 

Execution is terminated if this amount is not more than 20% greater than 
the amount required. READMX is called to read the MIXCLASS card and 
passes back the array of percentages. 

After the END card is read, checks are made to determine whether user 
requests are valid and conplete. 
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SMORDR'3 


Subroutines called by SMORDRt MMTAPE 

CTZMRD 

BRPRNT 


REAI»1X FVAL 

CTLPRM TSPACE 

IVAL 


3. Input Description 


Control cards are read by calls to CTLWRD. If a disk is used, it is first 
checked to make certain that the file exists. 

4. Output Description 

"Options chosen" messages are typed in addition to requests for more 
information. 

5. Supplemental Information 


Not applicable. 


6. Flowchart 



LARS Program Abstract 


MODOLB IDEWTIPICATION 


Module Mamet readicc 

Purpose t Inf gprets the MixciAi 

Sye tem/Langoage t cms/fortran 


Function Name: 


Date: 



PURDUE UNIVERSITY 
Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 
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READMX-2 


1« Module Dsaqe 
REMMX 

CALL READHX (LCARD, COL) 

Input Arguments: 

LCAAD - L*l, Card image of the MIXCLASS card being interpreted. 

COL - 1*4, t)ie column number preceding the first name on the 

MIXCLASS card (i.e., the first nonblank after the keyword). 

Output Arguments ! (passed in SMOOOM) 

NCLNAH- 1*4, NCLNAM is a two dimensional vector that is filled 
with the new class names from the MIXCLASS card. 

NEWCLS- 1*4, the nundaer of new classes. NEWCLS is incremented 
as each new name is read. 

MIXDA1V 1*4, MIXDAT is filled with the percentage ranges 

specified on the MIXCLASS card. MIXDAT (1,1) is the lower 
boundary of the percentage range and MIXDAT (1,2) 
is the higher boundary. 

2. Internal Description 

READ MX can be reduced to several functions — stripping off the name, 
scanning and modifying the user supplied ranges, and the reading and 
loading of the ranges into MIXDAT. 

Ihe first function is performed by using LOCATE to find the left parenthesis 
of the class being interpreted and then by a call to BCDFIL. Ihe data for 
this particular name is then scanned and the positions of all hyphens are 
noted by set«J.aq a corresponding flag in a flag vector. After a hyphen 
is read it is changed to a comma so that once the data is finished being 
scanned it can be read by IVAL. Once IVAI has been called the flag vector 
is used to transfer the data cad into the array MIXDAT. This is done by 
placing the first data value into MIXDAT (1,1) (i.e. the lower boundary 
of the user defined percentage range) and then checking the corresponding 
flag in the flag vector. If the flag was not set, then this lower 
boundary is sinplv copied into the upper boundary, MIXDAT (1,2) If the 
flag was set then the next data value is treated as the upper boundary. This 
is continued until all values have been loaded. 

Subroutines called by READMX: LOCATE 

BCDFIL 

IVAL 


3. Input Description 

READMX checks for various syntax errors on the MIXCLASS card. It is 
therefore capable of requesting a new MIXCLASS card by using CTIMRD. 
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READMX-3 


4. Output Deecription 

RBAONX oan write a syntax error massage to the user, nils gives the user the 
approximate location of the error and asks for the card to be retyped. 

5. suapiesmntal Information 
Not applicable. 

6. Flowchart 


Not c^lioable 
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LARS Program Abstract 


MODULE IDENTIFICATION 

Nodule Wareet smooth Function Name t smootoresults 

Purpose t Cbmpute the dominant class in each cell and reclassify 
Sys tent/Language t cms/portram 

Authors John Cain Date: 4/18/80 

Latest Revisors Dates 


MODULE ABSTRACT 


SMOOTH takes "n” lines of a Classification Results File (where n is the 
nustber of lines in a cell ) , determines whi^ group or class makes up the 
largest weighted percentage of each cell, and reloads the cell with 
that class. User^defined priority classes remain tmchanged. 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 

west Lafayette, Indiana 47906 
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SMOOTO-2 


1. JloduXe U 0 Me 

8H001M (NOCOIA) 

Ii»ut ArquBantsi 

N0C0L8 • I*4» MOOOILS is the number of columns (i.e. the 
buffer length) that is to be segmented into 
individual cell widths. 

CELCOL “ 1*4# the nunber of columns in the user defined cell. 
CBLROH - 1*4# the nunber of rows in the user defined cell. 

2. Internal Descriptiwi 

Upon entry# SMOOTH uses NOOOLS and the user defined parameters CELCX)L and 
CEUUW to determine the total number of cells to be processed. Using this 
information the main program loop is entered. The main loop begins by 
tallying the nunber of samples of each group in the cell. If the 
"weight” option is in effect then each grot^'s tally is weighted. The 
total wei^t is then found bv summing all the weighted values. Each 
group weight is then dividea oy the total weight to obtain a cell percentage 
for its group. These percentages are conpared to user supplied group 
percentages from the MIXCLASS card. If the values fall within the 
MIXCIASS ranges, then all the samples of the cell are changed to this 
MIXCLASS group and rewritten into the main buffer. Any user defined 
priority sanples (from the PRZORITy card) will not be altered.. If none 
of the MIXCLASS ranges correspond to the calculated percentages, then the 
nonpriority group with the greatest percentage must be determined. For 
efficiency the groips not found in the cell are removed from the percentage 
array (conpressed) . The array is then sorted so that the group with the 
greatest percentage is first and the one with the smallest is last. The 
number of the first group found that is not a priority class is then 
loaded into a pointer array. Groups defined as priority are loaded 
(overlayed) into their corresponding position in this array. The cell 
is then reloaded by using this pointer array for indexing. This loop is 
repeated for each CELCOL group of columns across the input Classification 
Results lines. When all cells have been modified, control is returned 
to 8M0INT. 

3 . Input Description 
Not applicable. 

4 . Output Description 
Not appliced>le. 

5. Supplemental Information 

The buffer of sacples modified by SMOOTH is passed through GLOCON in ARRAY. 
It's dimensions are CELROH by NOCOLS. 

6 . Flowchart 


Not applicable. 


• 
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LAKSYS COlSiTKOL CAKOS 
CHANbfcOiTeCTlON 


SSSoCCOL*!! ^AHAMEVt’R 


function 


default 


♦ •CHAN6EDETECT (NONE) 




CNONE) 


♦ BASERESULTS 


TAPeCXAA) 
FILE 
DISK 


IFF) 


0^. results mom (NONE) 

^^CAT|p 8 n tape XXA. 


(SEE control CARU) 
( dictionary ) 


♦ compareresults 


TAPE (XAA) 

FILE(FF) 

DISK 


location of RESULTS FROM 
SECOND DATE. 

LOCATED ON TAPE XXX. 

FILE FF. 

USE RESyLTS,PLACEU CN UISK 
IN CURRENT TERMINAL SESSION, 


(NONE) 


NEmRESULTS 


TAPE(AXx) 
FlLf (FF) 
INIT 

DISK 


LOCATION OF^CMANGE (OUTPUT) 
RESULTS FILE. 
ikRITE on tape XXX. 

FILE FF. 

inItialue Tape and write 


itialue tape 

RESULlS in FTl 
place results On 


E 1 
DISK 


(NONE) 


♦ HLOCk 


RUN(XXXXXXXX) 

LINES(x.Y«Z) 

COL(XtY*Z) 


RUN NUMrttM IS XXXXXXXA- 

niSPLAY lines X to y 

wITn line INTE*<YAL Z 
DISPLAY COLUMNS X TO Y 
WITH COLUMN interval Z — 


-I 


•I 

I 

1 (NONE) 


♦ DATA 


•START OF DATA DECK* 


define Change classes of 

ON A (CLASS* CAPUt AND 
FROM THE first 
WHICH CLASSES Fi<UM The 
•COMP* CARO» 

TO This class. 


interest bY NAMINo TH£ CLASS 


INOICATINO «iHlCH CLASSES 
CLASSIFICATION (ON *bASE* C 
The SECOND CLASSIF" ‘ 


SECOND 


ARD) AND 


_ SSIFICATION (ON 

permitted FOR A POINT TO BELONG 


class NAME! 
BASE Nit N2« 
COMP Ml* M2* 
CLASS NAME? 


Efc 


N3* 
M J* 


(WHERE Ml, M2*... AND Nl*N2*... ARE Cl ASS Oh 

(POOL NUMbtHS FNOm classifypoints. 


♦ END 


(NONE) 


END OF FUNCTION. 


(NONE) 
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Z<ARS Program Abstract 1058 


WODOLfe IDENTiriCATIOW 

Nodule Waiwt chasw Function Wamo t qiangedetectiom 


Purpooui Supar^sor for the CriANGE DETTCTIOn function 



PUROUe UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 
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CHASUP-2 


1 « Module .usage 

CAUi CHASUP 

There are no argiiments to CHASUP* It is called from LARSMN when the 
CHANGEDETECTXON function is requested. Control returns to LARSMN upon 
completion of the function. 

2. Internal Description 

OlASUP receives control from LARSMN to perform the CHANGEOETECTION processing. 
OiASUP calls CHAROR to read and interpret the control cards. Upon return 

from CHARDR* CHASUP calls CHANGE to finish the processing. Subroutines 
called by CHASUP: CHARDR 

CHANGE 

3. Input Description 
Not applicable. 

4. Output Description 

TWO informational messages, CHANGEDETECTION FUNCTION REQUESTED and 
CHANGE FUNCTION COMPLETED, are written at the typewriter. 

5. Supplemental Infor.-nation 
Not applicable 

6 . Flowchart 


Not applicable 


91 


LARS Program Abstract 1059 


MODULE IPENTIPICATIOW 
Nodule Name: chaoom 

Purpose : 


Function Name : changbdetection 


Block data 


System/Language i CMS/Fortran 
Author: John Cain 


Date: 6/1/79 


Latest Revisor: 


Date: 


MODULE ABSTRACT 


This is the BLOCK DATA subroutine for the CHANGEDETECTION common block CHACOM 


PURDUE UNIVERSITY 

Laboratory for Applications of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 




CHARDR interprets all function control cards tor CHANGEUETECTION. ChecKs 
are made for complete and valid specifications emd the proper input- 
output devices are attached. 


j 

i 


•» 

I 


PURDUE UNIVERSITY 



Laboratory for Applicatioas of Remote Sensing 
1220 Potter Drive 
West Lafayette, Indiana 47906 
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CHARDR«2 


1. ilo4wto U»aqe 
CALL CHARDR(Z,NAMB) 

OtttPttt Arquaantst 

2>L06ICAL*1 ea^ element initialized to .false. 

2(i,l» “ if class i from the base classification is part of 
user defined class j. 

Z(m«2*n)«.TRUB. - if class m from the 2nd (conpared) classification is 
part of user-defined class n. 

NANB-IM “ contains the names of the user-defined classes 

Listed below are the actions taken when the following control cards are read. 

BASERESULT8 TAPE - the variable TAPEl (TAPr.2) is set to the given 
(OOHPARERESULTS) tape nun^er. 

FILE - the variable FZLEl (FILE2) is set to the given file 
number. 

DISK - DZSKFG is checked to be sure that the DISK option 
is not already in effect, the tape and file numbers are 
cheeked to be sure that both the DISK option and TA«E option 
are not being used simultaneously. if they are, then an 
error message will be printed and the DISK will be used. 
RESLTl (RESLT2) is set equal to CXASSR. 

NCWRESULTS TAPE - the variable TAPE3 is set equal to the given tape 

number. 

FILE - the variable FILES is set equal to the given file 
number. 

ZNIT - the variable INITFG is set equal to 1. 

DISK - the same checks are made as above in addition to a 
dieck to see whether the INIT and DISK option were used 
simultaneously. DZSKFG is set equal to one and RESLT3 is 
set equal to CLASSR. 

BLOCK RUN — the variable RUNNUM is set equal to the given run 

number. 

LINE - STALIN is set equal to the first entry (the 
starting line of the area to be investigated) . LASLIN 
(last line) is set equal to rhe second entry and finally 
LINZNT (line interval) is set equal to the last entry. 

OOL - same as above where the variables are: STACOL - 
first entry, LASCOL - second entry and COLINT - final entry. 

DATA A check is made for the presence and validity of all 

information. 

CLASS name Ihe name given is stored in the array NAME. 


CHAROR'3 


BASE 

OOMP Ml, M2... 


9A 


Usinq the given dess nunbere the appropriate locations in 


the 2(64^2,64) array are set .TRUE. ^ 

(i.e. if these are the BASE and COMP cards for the j 
user-defined CLASS, then the following assignments are 
made for array Zi 

Z(Nl.l.j) • .TRUE. 

ZOa.l.j) • .TRUE, 
and r 


Z(M1.2,j) » .TRUE. 
Z(M2.2,j) • .TRUE. 


2. Internal Description 

OIARDR uses the standard card reader logic in using CTLMRD, CTLPRM and IVAL 
to read and interpret control cards. 

QIARDR begins by initializing all flags and arrays that are used to convey 
control card information. It then goes into a loop of reading and interpreting 
the input specifications and the BXOCK card, vmen the DATA card is read 
CHARDR checks for the presence of all information and its validity. Another 
loop is entered and the CLASS cards and their corresponding BASE and COMP 
cards are read. %e class nundbers from the BASE and COMP cards are used to 
set appropriate values in the Z array to a logical .TRUE. 

Z(i,l, j)*.TRUE. if the Ith class from the base results file is part of user- 
defined class j. 

Z(k,2«m)«.TRUE. if class k from the conpared results file is part of user- 
defined class m 


Ihis loop is exited when an END card is read. Once this card is read. CHARDR 
calls CHTAPE to mount the specified tapes. If a disk was specified as an 
input device. CHARDR first checks to be certain both a tape and disk were not 
specified for a single input. It then reads from the results file to be sure it 
exists on the disk. If a disk was specified as an output device, checks are 
made to be sure there is sufficient space for the output results. TSPACE maJ^es 
a search for a larger disk if necessary. CHARDR finally returns control back to 


CHASUP. Subroutines called by CHARDR: 



CTLWRD 

CTLPRM 

TSPACE 

BCDFIL 

CHTAPE 

RTMAIN 

IVAL 

ERPRNT 



3. Input Description 

Function control cards for CHANGEDETECTI ON are read via CTLWRD. 

4. Output Description 

control card error messages are written via ERPRNT. 

5. Supplemental Information 
Not applicable. 

6. Flowchart 


Not applicable. 
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LARS Program Abstract 


WOPOLB tPEWyiPICATION 

Modioie Wamoi chanse Function Mama : changedetectiow 

Porposot Cowoaras t«ID classification regultg filoa unrt 

SyPtoiK/tafngttage * CMS/Portran 

Atfthor : 3dhn Cain D ate : fi/3/79 

Latest Revisor: ^Date: 


MODULE ABSTRACT 


CHANGE is the main subroutine for CHANGEDETECTION. It reads from two input 
tiqpes (or one disk and one taqpe ) , calls COMPAR, then outputs the data in 
standard LARSYS classification results file format to tape or disk. 
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CHANGE- 2 


1 » Modttle Usage 
CHANGE 

CALL OIANGE(Z,NAME) 

Input Arqueantst 

Z-Logical*l - Z(i,l, j)«.TRUE. if class i from the base classification is 
part of user-defined class j. 

Z(m,2,n)*.TRUE. if class m from the 2nd (compared) 
classification is part of user-dufincd class n. 

NAME - X*4 - contains the names of the user-defined classes. 

Output Arguments; 

Mot applicable. 

2. Internal Description 

CHANGE first reads the file numbers from the input tapes, and the tape 
nunfi>ers passed through the common blocl^ and creates a code that takes the 
place of the CLASSIFICATION STUDY number. The code format is the base tape and 
file numbers followed by the comoared tape and file numbers. CHANGE then reads 
the area identification record (record type 5) from both input sources and 
checks to see whether they are valid for the given BLOCK CARD; if not, 
appropriate error messages are printed and the function is terminated. Record 
types 1-5 are written to the output tape (DISK) . The inputs are positioned 
to the correct line number and shifted to the correct column number. CHANGE 
then calls CO.MPAR to determine which class each point belongs to and this 
information is used to create file type 5. Finally record types 7 and 8 are 
written and control is returned to CHASUP. If the output device is a tape, then 
a final record type 1 and END OF FILE Mark are written before returning to 
the supervisor. Subroutines called by CHANGE: 

COMPAR 

RTMAIN 

TAPOP 


3. Input Description 

Record types 1, 5, 6 of the LARSYS classification results files are read 
from the two input devices, RESLTl and RESLT2. One of these may be a disk 
(DSRN C..ASSR) . Tape drives 181 (CPYOUT) and 182 (CO.*!PTP defined in CHACOM) 
are used as inputs. 

4. Output Description 

The output device RESLT3 initially has a DSRN of MAPTAP. If a disk is used 
(only if one is not used for input) , the DSRN is changed to CLASSR. Tape 
drive 180 is used for output to facilitate the run of prtntkesulTS on the 
output data immediately after the CHANGEDETECTI ON run. The output is a 
classification results file in standard LARSYS format. 
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5« SttPPlanintal InfonaatiTO 
Mot applioablo. 

6. Plcwonart 
Not applioablo* 


CHANGB-3 
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LARS Program Abstract i062 


MODOLB IDENTIFICATION 

Modulo Naraet compar Function Mama t chance 

Purpose : Compats 2 lines of classlflcatlwi results 

System/Language t cms/fortr^ 

Author : Susan Schwlngendorf D ate : 3/28/79 

Latest Revisor t Date: 


MODULE ABSTRACT 


COMPAR compares two lines of classification results (presumably from 
two different classifications which are registered to each other) against 
user defined change classes in a logical array* and writes the output class 
nund>er in the output vector. 
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COMPAR-2 


1. tipdttl* ItoAM 

CALL COMPAR (NGOLS, NCLA88.Z, BUFFI, BUFF2. BUFF3) 

IttPttt Argument ■ ; 

NOOLS “ INTEGERA4, the nuniber of columns of classified data* 

NCLASS - 1NTCGER*4, the number of classes defined by the user 
in array Z. 

Z *> LOGICAL*! (64, 2, 64) array containing user defined 
change classes. (Initialised to .FALSE.) 

z(‘:2;k) - '*“• 

I from classification 1 (BUFFI) and in class 
J on classification 2 (BUFF2) should be assigned 
to class K in BUFF3. 

BUFFI - LOGICAL*! (2*NC0LS + 4) vector containing classified 
data from first classification. First full word is 
line number. Then the second byte of each halfword 
contains the next class number 

BUFF2 - LOGICAL*! (2*NC0LS -i A) vector containing classified 
data from the second classification. First full word 
(4 bytes) is the line number. Then the second byte 
of each halfword contains the next class number. 

Output Arguments ; 

BUFF3 - LOGICAL*! (2*NCOLS + 4) vector of change classes for 
this line. The first full word contains the line 
number. Then the second byte of each halfword 
contains the assigned change class number. 


2. Internal Description 

The line number is written in the first word of BUFF3. The next class 
number is then extracted from BUFFI and BUFF2 and assigned to integer variables 
CLASSl and CLASS2. A loop through the logical array Z determines which output 
class to assign this point to. If Z (C1.ASS1, 1, J) and Z (CLASS2,2,J) are true, 
then the point is assigned to class J. If it belongs to none of the defined 
output classes, then it is assigned a class number NCLASS+1. The output class 
numbers are written in BUFF3. 


3. Input Description 
Not applicable 

4. Output Description 
Not applicable 

5 . Supplemen t al Information 
Not Applicable 
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CONPAR-3 


Flowchart 
Not Applicable 
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LARS Program Abstract if>fe3 


MODULB IDBNTiriCATIOW 
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Author * E.M. Rodd ^Da te : 9/5/ 72 
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MODOLE ABSTRACT 

CHTAPE motints and positions the results tape (or a tape to be used as output 
for copying results files) . 
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CHTAPE-2 


1. Hodul« 0»Aq« 

CHTAPE 

CALL CHTAPE (RQTAPE,RQFILE,MODE»UNIT) 

Input Argunwintst 

RQTAPE - 1*4, Tape number of requested tape. A 

tape nuinber of 0 is a request for a scratch 
tape. 

RQPILE •* 1*4, File number of requested file. If 
RQFILE is « 0, then the tape will be 
initialised by writing a record type 1 
on the results tape with filetype • 0. 

MODE - 1*4, Flag indicating usage of CHTAPE. 

MODE ■ >1 indicates CHTAPE has been called to 
mount and position a tape to be used for 
copying results files onto. Mode « 0 
indicates that a results tape is being moxmted 
for reading a results file. In this case, the 
tape is mounted ring out. Also, if MODE ^ 0, 
RQFILE « 0 is invalid and will cause an error 
when an attempt is made to write on the tape. 
MODE * 1 indicates a tape is being mounted 
for writing a new results file (or continuing 
a suspended classification) . The difference 
between MODE « -1 and MODE ■ +1 is the DSRN 
\ised for the tape. For MODE « >1 DSRN is 
CPYOUT and for MODE » tl, DSRN is MAPTAP. 

(DSRN is MAPTAP for MODE > 0) . 

UNIT - 1*4, DSRN of tape being mounted. 

Output Arguments ; 

RQTAPE - 1*4, When MODE « 0, set to -1 if requested 
tape file was full and user decided to use 
disk for results. Otherwise, remains 
unchanged. 

RQFILE - 1*4, When MODE « 1, set to -1 if requested 
tape file was full and user decided to use 
disk for results. Otherwise, sends back 
current file position of tape. 

CHTAPE checks the validity of the tape by reading the record type 
1 from the tape and verifying the tape and file number as well 
as checking for the correct type of file. Any 
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CHTAPE-3 


attempt; to overwrite an existing file causes CHTAPE to ask 
I the user (via the typewriter) if he wishes to overwrite the 
lfile« respeeify a new results card, or terminate the function. 

Note, however, that if a request has been made to initialize a 
tape, AO checking is performed on previous contents. 

2. Internal Description 

See Output Description. Subroutines called by CHTAPE: 

TAPOP RIN6IN IVAL 

MOUNT CTLWRO ERPRNT 

CPPUNC CTLPRM RTMAIN 

3. Input Description 

The record type 1 of the results tape is read for each file 
up to and including the file needed. That is, if file 4 xs 
requested the record type 1 is read from files 1-4. 

4 . Output Description 

The following information messages are issued under the 
circumstances listed. The term filetype means the filetype 
code from record type 1 of a results file (the program uses 
variable CHECK for this number) . 

I0C42 is typed when a tape has been mounted and before 
CHTAPE positions it. This message is not typed 
when the tape is being initialized or when the correct 
tape number was already mounted. 

10043 is typed when MODE » and filetype of the requested 
file *0. 

10044 is typed when MODE » 4l and filetype of the requested 
file s 1 and the restart flag fromGLOCOM (RESTRT) is 
not * 1. 

10045 is typed when the tape is correctly positioned. This 
is not typed when initializing a tape. 

After 10043 and 10044, the user is asked whether he wishes 
to overwrite the file, respecify a new results card with a 
I new tape and/or file or disk option, or terminate the function. 

10100 is typed to allow entry of the new results card. 

This occurs when the user requests to respecify the 
results card. 

10101 is typed to confirm usage of disk for results and 
occurs whenever disk is specified on the results 
card. 
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CHTAPE-4 


The folloiiring error messages are typed under the conditions 
listed. 

E361 is %nritten when the tape Is being filed forward and 
a file is encountered with filetype other than zero 
before the requested file is reached and MODE « 0. 

E362 is written when the circiunstance for E361 occurs 

and MODE « 1. It is also written when MODE “ 1 and 
the filetype of the file requested is s >i. 

E363 is written if the RESTRT flag is « 1 and the filetype 
of the requested file is not = 1. 

E364 is written when MODE “ 1 and the filetype of the 
fiel requested « 1. 

E365 is written when an EOF is read on the results file. 
This should never occur with valid results files. 

For message texts refer to the User's Manual. 

5 . Supplemental Information 

This section deals with the handling of tapes by CHTAPE 
Input ; 

If a tape is mounted on the device and it is the incorrect 
tape number (as noted from the appropriate status words in 
GLOCOM ) t TOPRU is called to unload the tape before the 
correct tape is mounted. If the correct tape is mounted, 
CHTAPE will check for the ring in if MODE * tl. If the ring 
is not in, the tape is unloaded and MOUNT is called to 
mount the tape with the ring in. If the correct tape is 
mounted, CHTAPE assumes that the file number (as recorded 
in GLOCOM) is correct and moves the tape backwards or 
forwards to find the requested file. 

Output : 

The tape is moxmted with ring in for MODE « ^1 and with ring 
out for MODE > 0. 

The tape is left positioned at the beginning of the requested 
file. When the tape is initialized a TOPRW is used to do 
this. 

6 . Flowchart 


Not Applicable 
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Exhibit 1 


FRIS III Timeline Chart 


'^«**'* Technoudgy Transfer 


Activity: 


A. Training 

1. Short Courses 

2. VbRKSHOPS 

3. Rioto- Interpretation 

Sort Course 

B. Consultation 

C. Documentation 

1, LARS User Documentation 

2. NCC User Dxufentation 

D. Terminal Operations 


Calendar Year 


1979 


4/1 - 
6/30 


7/1 - 
9/30 


10/1 - 
12/3 


SL 




1980 


1/1 - 

-JZ3L 


4/1 - 
6/30 




.S7 


^ 


Ji SM 


V 

▼ 

SL 


planned start of activity 
actual start of activity 
duration of activity 
progress toward activity completion 
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Exhibit 2 


FRIS III Timeline Chart 







FRIS LARSYS Software Documentaticm Task 

























856C 











Jun I Jul I Aug 
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ISC r*f* tisac <••? ia« ai> «msa.isc 































Exhibit 3 


FftlS III Timeline Chart 


Calendar Year 


Task: PREPRKESSING TRANSFER 


Activity: 

A. Hanning 

B. Program Hfinement 

L Landsat 3 Reformahing 

2. Geotctric Correcticn 

3. Image Registivvtion 

C. Program Transfer 

D. Consultation & Debugging 

E. DoourenATiON 

F. Test & Bmluation 

G. SiPPORT Activities 

1. Lancsat 3 Data Evaluation 

2. FRIS N’ Coordinates Definition 

3. teFORMATTiNG Operations 

Procedures 


4/1 - 7/1 - 10/1 - 

6/30 9/30 12/3 



planned start of activity 
actual start of activity 
duration of activity 
progress toward activity completion 


V 

▼ 







3 


Exhibit J! 

FRIS III Timeline Chert 


Task: fWMfiErenr 


Activity: 

A. f^PORTING 

1. Informal Monthly Status 

2. Monihly Fiscal 

5. Quarterly Progress 
A. S^i-Annual I^ens 

B. Information Dissemination 

C. Cost Evaluation 

D. SPECiM. Projects 

L Cu^iFicATioN Accuracy 
Evaluation 

2. Ratio Evaluations 

3. l^iABB Application Test 


1979 

4/1 - 
6/30 

7/1 - 
9/30 

10/1 - 
12/3 


planned start of activity 
actual start of activity 
duration of activity 
progress toward activity completion 










